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Summary. We review some theoretical aspects connected with the interplay of 
strong electron correlations and vibron phenomena in negative fuUerene ions and 
in solid fuUeride conductors, superconductors and insulators. The first part cov- 
ers molecular ions, their intra-molecular Jahn- Teller effect, Coulomb (Hund's rule) 
exchange interactions, molecular vibrons and multiplet splittings. The second part 
addresses electron propagation in molecular fullerides, with special emphasis given to 
trivalent cases such as KsCeo and NH3 K3C60, where metallic, superconducting and 
Mott insulating phases are at play. Dynamical mean field theory approaches to a sim- 
plified Hamiltonian for this system are discussed in the light of some of the observed 
phenomenology. It is argued in particular that NH3 K3C60 is a Mott- Jahn- Teller in- 
sulator, which under pressure turns into a strongly correlated superconductor, thus 
sharing some important elements with the high-Tc cuprates. 



6.1 Introduction 

The wealth of experimental data on the alkali fullerides collected through the 1990s 
and early 2000s has raised several exciting theoretical issues which, in turn, have at- 
tracted a large research effort. Different phases are realized when temperature, alkali 
concentration, lattice parameter and lattice structure are varied. Metallic, insulat- 
ing, and superconducting phases are obtained in solid compounds characterized by 
the partial occupancy - by an average n electrons per Ceo with n ranging between 
and 6 - of the narrow electronic band originating from the LUMO molecular 
orbital. When this band is completely empty (pure Ceo) or completely filled (KeCeo) 
the solid is a band insulator (although the lattice structure in the two solids are dif- 
ferent). Band-structure calculations of compounds characterized by an incomplete 
filling of this band consistently predict metallicity [1, 2, 3]. However, experimen- 
tally some of these compounds (e.g. K4Ceo, Rb4Ceo, NH3K3Ceo) are found to be 
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insulators [4, 5, 6, 7, 8, 9, 10]. Among the metallic compounds, many with n = 3 

turn superconducting at low temperature, with transition temperatures as high as 
30 — 40 K, while several others remain metallic at all experimentally investigated 
temperatures. 

For even n, as in K4C60, electron-electron correlations and the JT coupling sta- 
bilize a correlated insulating state of a lattice of evenly-charged Ceo molecules. In 
this type of insulator, fluctuations about (ui) = n are suppressed, and a gap opens 
in the electronic spectrum. This state is non-magnetic, very much like a regular 
band insulator. However, the mostly intrarmolecular electron correlations responsi- 
ble for band narrowing and gap opening axe largely coulombic and vibronic in origin. 
We have suggested that the (body-centered tetragonal [11]) structure of K4C60 and 
Rb4C6() may be a realization of this state where electronic and vibronic interactions 
play an important role. Similarly, the insulating state observed in strongly corre- 
lated NH3K3C60 (and (NH3)6Li3C6o [12]) as well as the insulator-superconductor 
transition obtained under pressure are likely to have the same origin. 

A satisfactory understanding of how similar physical parameters lead to very dif- 
ferent ground states (insulator/metal/superconductor) has long eluded the research 
community. A novel picture relating the transitions between metallic, insulating, 
and superconducting phases to the interplay of strong electron-electron correlation 
with electron-phonon coupling in the LUMO band has emerged in recent years. The 
scenario is now in our view better understood, even though many quantitative de- 
tails still escape the full grasp of theory. The present Chapter reviews and illustrates 
this picture. 

6.2 Ceo molecules and molecular ions 

We review for a start some basic aspects of the physics of Ceo , concerning electron- 
electron correlation and electron-vibration coupling in Ceo molecular ions and ionic 
fullerides. For more complete reviews we suggest consulting Refs. [13, 14, 15, 16]. 

In Ceo all the 60 carbon atoms of fuUerene are arranged as a regular, icosahedral, 
roughly spherical cage, of ~ 0.7 nm in diameter [17]. This cluster may be thought as 
a piece of graphenc sheet, wrapped up to a spherical shape. The regular hexagonal 
structure of graphite is distorted, with 12 five-membered rings intercalating 20 six- 
membered ones, thus introducing the Eulerian An solid angle necessary to yield a 
closed surface. The order is such that all carbon atoms remain equivalent, each 
sitting at the corner of one 5-member ring and two 6-member ones. 30 chemical 
bonds, the so-called double bonds, arc shared by two hexagonal rings only; the 
remaining 60 bonds, shared by pentagonal and hexagonal rings, arc about 5% longer 
[18]. The cr-bonding sp^ graphite orbitals constitute the backbone of the molecule. 
Spherical curvature alters their character to an average "sp^'^*" [19]. In energy, the 
(T-bonding orbitals range from several eV to a few tens of eV below the vacuum, the 
antibonding states lying -1-10 eV and higher above vacuum zero [20]. The chemically 
active electronic states are those derived from the "pz" carbon orbitals (actually 
of hybrid "s^'^^p" nature [19]) that are directed radially, supporting a half filled 
TT-electron system. 

Although theoretically metastable compared to graphite and diamond. Ceo is 
nevertheless an empirically very stable and long lived allotliropic form of carbon 
both as an isolated molecule and as an fee, or Pa3, solid. The reason for stability is 
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the substantial amount of electronic energy gained in the delocalized MO's, making 
in some sense Ceo the spherical counterpart of the aromatic benzene ring. The 
molecular orbitals at the origin of the conduction bands of solid fuUoridos are all of 
TT-bonding/antibonding nature. 

The overall molecular symmetry group is the icosahedral group Th, the largest 
point group in 3D (except for axial groups). Symmetry implies a large degeneracy of 
the group's irreducible representations [21]: Ag/u (1-dimensional), Tig/^, T-zg/u (3- 
dimensional), Gg/^ (4-dimensional) , and i^g/„ (5- dimensional). Large degeneracies 
are accordingly very common to all electronic, vibrational,vibromc molecular states. 
This makes Ceo a rich playground of novel vibronic structures, where Berry phases 
and entanglement plays a fascinating role [22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32]. 

Icosahedral symmetry is of course restricted to the ideal situation of a molecule or 
molecular ion in vacuum. In compounds and/or in the solid state it will be affected 
by crystal fields of lower symmetry. Even in vacuum, isotope substitution of one 
carbon atom is enough to reduce Th to simple bilateral reflection Cv Due to 1.10% 
isotope abundance of ^"'C in natural carbon, only about 50% of Ceo is pure ^'^Ceo. 
However, isotope substitution induces only small splittings (about 1%) of the vibron 
frequencies, and can be safely neglected for many purposes. For a discussion of the 
intricacies of isotope shifts in solid-state properties of the fuUerides, see Ref. [16]. 

As mentioned, immersion of the icosahedral molecule in a solid-state environ- 
ment reduces its symmetry to that of the local crystal field. For example, solid Ceo 
has electronic bands and optical plionons compatible with the local cubic field in- 
duced by the fee lattice. In most (but not all) solid state compounds, intermolecular 
interactions are relatively loose, each molecule retaining its close structure. Accord- 
ingly, splittings of the molecular vibrations (now optical phonons in the solid) due 
to reduced symmetry are small and hardly observed at all [33]. 

6.2.1 Molecular electronic states 

Electrons in the tt orbitals of Ceo represent the chemically relevant region of the 
spectrum. These orbitals provide the basic one-electron picture, neglecting first e-e 
correlation effects, which we shall address in Sec. 6.2.4. 

Many approaches have been taken to the electronic structure of Ceo, from simple 
Huckel tight-binding with one orbital per atom (already yielding the correct order 
of molecular energies and gaps) [20, 34, 19], to more extended bases [20, 35] to mi- 
croscopic DFT-LDA (density functional theory in the local density approximation) 
calculations on localized [20, 36] and extended bases [1, 37]. 

These approaches provide an increasing degree of quantitative accuracy in the 
description of the electronic spectrum. Simple intuitive models, such as the particle- 
on-a-sphere model [38], are often sufficient for a qualitative and synthetic under- 
standing of the electronic structure. In this approximation the idea is treating the tt 
electrons as though moving on a spherical shell of radius R, mimicking the attrac- 
tive potential generated by the carbon ions. The precise localization of these ions 
produces then a weak icosahedral perturbation. As long as this icosahedral pertur- 
bation can be neglected, the single-electron angular wavefunctions are atomic-like 
spherical harmonics Ylm, with energy 
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Fig. 6.1. A particlc-on-a-sphcre schematic representation of tlic electronic levels 
of Ceo- The HOMO and LUMO levels, originating from the L = 5 orbital, are 
indicated. (Adapted from Ref. [38].) 



where me is the electron mass. 50 out of 60 Pz electrons of the neutral molecule 
fill completely the MO up to i = 4. The lowest L = 0, 1,2 orbitals coincide with 
icosahedral states labelled ag, tiu, hg respectively. All higher L values are split into 
icosahcdral representations by the icosahedral field generated by the cage. After 
filling all states including L = 4, 10 electrons are left in the L = 5 shell, which 
is therefore only partly filled. As it happens, the icosahedral splitting (L = 5 — » 
hu+tiu+t2u) of this 11-fold degenerate orbital generates a closed-shell configuration, 
as shown in Fig. 6.1. In accord with microscopic calculations and with experiment, 
the completely-filled HOMO has hu symmetry, and the LUMO is tiu- The HOMO- 
LUMO gap is therefore caused by the icosahedral perturbation in the L = 5 shell, 
and is experimentally ~ 1 eV for molecules in vacuum [39]. A t^g LUMO-l-1 state, 
originated from the L = 6 shell, is found approximately 1 eV above the tiu LUMO. 

The electron affinity of Ceo is large (2.69 eV) [40, 41] and experimental evidence 
has been found that CgQ [42] and even Cgg [43] are stable ions in vacuum. In solution, 
a wider spectrum of ionization states has been demonstrated electrochemically, up 
to Cg^ [44, 45, 46, 47]. As an adsorbate on a metal surface, the electronegative 
Ceo molecule naturally picks up electrons [1, 48], and evidence has been provided 
of charge transfer as large as n=6 [49]. In the solid state, compounds have been 
synthesized, covering a wide range of charge transfers, from n—1, as in TDAE-Ceo 
[50, 51] or RbiCeo [4], n=3, as in KgCeo or RbsCeo [52], n=4 as in K4C60 [53], 
n=6 as in RbeCeo [54], or even higher as in Lii2C6o [55, 56]. More recently, also 
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Fig. 6.2. Spherical resolution of the vibrational spectrum of Ceo- The modes axe 
organized according to the three series, indicated with a;i/2/3- The spherical parent 
L is indicated before the icosahedral label. (Prom Ref. [65]). 



positive Ceo ions were produced in gas phase [57], in liquid solution [58, 59], in Ar 
matrix [60], and in the solid state [61, 62, 63, 64]. Since the LUMO can hold up to 6 

electrons, the negatively charged ions CgJ" with 1 < n < 5 are open shells. Likewise, 
the cations Cgg" with 1 < n < 9 are also open shells due to the fivefold degeneracy 
of the HOMO orbital. 



6.2.2 Moleculcir vibrations 

Ceo has 174 vibrational degrees of freedom, but thanks to symmetry-induced de- 
generacy and selection rules, the vibrational spectra show relatively few peaks, with 
clear identification. In particular, only the 4 Ti„ dipolar modes axe infraxed active, 
and only the 2 Ag and 8 Hg modes are Raman active. Neutron experiments [66] 
are sensitive to all modes, including the silent ones, but with low resolution. Cal- 
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culations are therefore crucial, in order to get a global picture of the vibrational 

spectrum [65]. 

Like for the electronic states, many routes to the calculation of vibrational cigen- 
frequencies and normal modes have boon pursued, both based on force field fits to 
the experimentally accessible data [66, 67, 68], or based on ab initio calculations of 
the molecular structure [35, 37, 69, 70]. The agreement among different calculations 
is somewhat worse than for the electronic levels: typical discrepancies on the order 
of several meV are well above nowadays' experimental resolution. 

A qualitative understanding of the vibrational structure of Ceo is provided by 
the the analogy of the Ceo cage with a hollow elastic sphere proposed in Ref. [65] and 
illustrated in Fig. 6.2. The cigonmodes of a homogeneous spherical membrane with a 
stretching and a bending rigidity are collected in three classes: mainly radial, mainly 
tangential and purely tangential. The first class contains levels of parity (—1)^, with 
L = 0, 1, 2, 3, the levels in the second class have the same parity, but start off at 
L = 1; in the third series the parity is reversed (—1)^"*"^, and they start at L = 1. 
The modes of rigid translation and rotation of the sphere are identified in the L = 1 
states at zero frequency in series 2 and 3. Of course, when the homogeneous sphere is 
replaced by the discrete 60-atoms molecule, the infinite set of spherically symmetric 
eigenmodes goes into a finite number of modes, now labeled by In representations. 
L = 0,1,2 states have icosahedral counterparts in Agj^, T-^^/g, Hg/^ respectively: 
the first gerade/ungerade label corresponds to series 1 and 2, the second labels the 
third series, which has inverted parity. States with L > 2 arc split, according to 
the rules given in Table III of Ref. [65]; for example a L = 3 vibration becomes 
T2u ® Gu, if its parity is odd. The explicit eigenmodes, computed for example by 
force-field methods, can be easily analyzed in terms of the spherical basis, to obtain 
their parentage in terms of hollow sphere modes [65]: the three series are readily 
identified in Fig. 6.2. Inter-mode mixing is present, but most vibrations hold a well- 
defined spherical parentage, so that the spherical picture remains generally valid 
and useful. 

6.2.3 Jahn- Teller coupling between electrons and vibrations: 
molecular vibronic states and energetics 

The electronic state of most Cg * ions is orbitally degenerate. According to the Jalm- 
Teller ( JT) theorem, the highly symmetric Th geometry is unstable toward symmetry 
reducing molecular distortions. Theoretically, such distortion can be evaluated by 
direct calculation of the relaxed geometry of the molecular ions, and are very well 
defined. However, direct experimental evidence of the JT distortions is rather poor, 
and limited to few compounds, mainly "discrete" salts, as discussed in detail in 
Ref. [15]. The reason for this is twofold: (i) the atomic displacents from the Ih 
positions of neutral Ceo are very small - of the order of few pm, and (ii) quantum 
tunneling and/or thermal hopping of the molecule between equivalent JT distortions 
generally restores, on average, the original Ih symmetry. 

Despite the small absolute atomic displacements, the JT physics involves serious 
energetics, since the phonon frequencies are high (owing to the large stiffness of the 
fullerene cage and the light carbon mass) and the JT couplings are large. At zero 
temperature, a significant fraction of the JT energy is associated to nonadiabatic 
effects related to the quantum motion of carbon nuclei, and cannot be obtained by 
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standard adiabatic calculations based on classical atomic positions. However, quan- 
titative insight in the coupled electron-phonon dynamics is obtained by means of 
the study of formally simple models. To the extent that these distortions are not 
too large, the distortion mechanism is well described within a model where the dis- 
tortion is expanded on the basis of the normal modes of vibration: these are treated 
as harmonic oscillators, and the coupling between the degenerate electron level and 
the JT-active distortions can be assumed to be linear in the phonon coordinate. 
These simplifying approximations remain valid only as long as the distortions are 
relatively small, as in Ceo ions, and could not be applied to ionic states of e.g. SiH4, 
whose molecular shape is radically modified upon charging [71]. 

The basic model Hamiltonian for the JT dynamics of an icosahedral molecule 
has the following standard structure [72]: 

H = Ho + He-. , (6.2) 

where Ho describes the free (uncoupled) electrons and vibrations and Hc-v the 
linear JT interaction. To begin with, consider a single vibrational mode of energy 
fiu! and symmetry label A: the Hamiltonian are written 

Ho = hiuY, (^M^ + ^) + (e - A«) cJn,aCm,a , (6.3) 

^e-v = gru,^ Yl H C'aX {bl + b^) ^ ^-m. o ■ (6.4) 

Here, the harmonic oscillator is represented by the boson operators 6jj: the indexes 
H span the degenerate representation A. The distortion coordinate appears in the 
coupling term (6.4) as Q/j, = 2~^^'^ (bj^ + 6,^^ . X = tin is the symmetry label of the 
LUMO electronic state, span by indexes rrii and generated by the fermion operators 
Cm a- is the dimensionless parameter measuring the electron-vibration coupling 
strength. C^^^ ^® ^ Clebsch-Gordan coefficient for the icosahedral-group [73], 
which recouples the fermion tensors of A symmetry with the boson tensor to an 
icosahedral scalar. In practice, this recoupling is possible only for modes of symmetry 
A — Ag, Hg. All other A's yield vanishing Clebsch-Gordan coefficients, which means 
that all modes of non-Ag / Hg symmetry are not linearly coupled to the tiu LUMO. 

It is interesting to note that the Ih Clebsch-Gordan coefficients for the coupling 
of two tiu representations to A = Ag/Hg happen to coincide numerically with 
the corresponding C;^^ Clebsch-Gordan coefficients for the full rotations group 
S0(3). The tiu label maps to angular momentum / — 1, while Ag maps to L = 0, 
and Hg maps to L = 2, in accord with the spherical picture of Figs. 6.1 and 6.2. 
This means physically that linear JT coupling of the tiu orbital does not distinguish 
between a soccer-ball molecule and a perfectly spherically shaped one: the tiu level 
behaves as an atomic-like p state perturbed by Ag/Hg distortions of monopolar 
and quadrupolar nature [72 , 74, 75]. The implications of this extra symmetry are 
striking: Hamiltonian (6.2) is actually a scalar with respect to the full rotation group 
S0(3), its eigenstates being representations of this larger group. S0(3) states are to 
be expected in the outcoming spectrum, which means that vibronic states belonging 
to different Jh representations collapse to degenerate SO (3) multiplets (labeled by 
some angular momentum L, rather than by Th labels), as on the left side of Fig. 6.1. 
In actual Ceo, the linear electron- vibration coupling term is just the leading term in 
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an infinite expansion, where only higher-order terms introduce the actual icosahedral 
symmetry, thus splitting the highly degenerate spherical vibronic states not unlike 
Fig. 6.1. Eventually, the effect of these higher-order terms is small, yielding very 
small splittings of the vibronic states: neglect of all these higher-order effects is a 
good approximation for CgQ~. 

Electron coupling to the A = Ag mode represents a trivial shifted oscillator, since 
the nondegcnerate mode does not split the electronic degeneracy, and the distortion 
is purely proportional to the total charge in the coupled level. Henceforth we shall 
mostly concentrate on coupling to degenerate A = Hg modes. 

There is no loss of generality in choosing the energy reference so that e = /i in 
Eq. (6.3). We are then left with a two-parameters Hamiltonian operator. The value 
hu! of the frequency of the harmonic oscillator sets the energy scale. Hq and -ffo-v 
are written with a common factor hui, which sets the energy scale of the model. The 
dimensionless parameter g tunes the intensity of the JT coupling, thus of the ten- 
dency of the system to distort. This is seen as follows: express each pair of vibrational 
operators 6^) as a dimensionless coordinate and conjugate momentum P^. 
In this notation, the adiabatic potential is obtained by ignoring the \ ^P^ phonon 
kinetic term in Ho and treating the operators as classical variables Qi^ [76]: 
one readily finds that Qo = g (all other = 0) is a minimum of the sum of the 
competing Q,j-liuear term (Hc-v) and Qfj-quadratic term in Ho. We see here that 
g measures the amount of distortion (in dimensionless oscillator coordinates) along 
the normal mode. For small g <S 1 (weak-coupling regime) distorsions are small. 
In this regime quantum fluctuations dominate: the correlated vibrational and elec- 
tronic dynamics involves all electronic states in the multiplet at the same time, in a 
profoundly nonadiabatic fashion. Despite this intricacy, as a small parameter can be 
identified, the JT problem can be dealt with by treating Hg-v by perturbation the- 
ory. For intermediate g ^ 1, the distortion is sizeable, but quantum kinetic energy 
has a relevant role in promoting tunneling among equivalent minima. The vibronic 
spectrum in this region shows nontrivial structures. Eventually, when the coupling 
becomes very large {g ^ 1), the JT distortion is so large that individual JT minima 
remain well separate so that tunneling is efficiently suppressed: the system freezes 
in one of the equivalent local minima. 

Here we should note a peculiarity of the so-called linear t (g) H JT problem 
Eq. (6.3)-(6.4) at hand, and precisely that the minimization of the lowest adia- 
batic potential sheet leads to a continuous manifold of equivalent minima (a trough) 
rather than isolated minima [22, 23, 72]. This is a consequence of the extra spherical 
symmetry of this special icosahedral linear JT system. However, at strong coupling 
higher-than-linear terms, neglected in Eq. (6.4) become relevant. They produce a 
warping of the JT trough, leading to isolated minima. The molecule distorts to a 
minimum in a set of static JT reduced-symmetry configurations separated by saddle 
points of the adiabatic potential surface. For the intermediate coupling case of Cqq, 
detailed Hartree-Fock calculations have shown, for example, that the total energy 
lowering in going from the X/j-symmetric molecular configuration to static JT distor- 
tions oi D^d, Dsd, D2h symmetries, are in fact identical to within 1% [77]. Therefore, 
the linear JT coupling approximation Eq. (6.4) is very well justified for Ceo ions. 

Analytical expression for the t ® H JT energies can be obtained in the simple 
weak-coupling limit, by means of perturbation theory [23]. The main outcome of the 
weak-coupling regime is that the ground-state energy gain E"'^ in this JT system 
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Fig. 6.3. Ground-state energy lowering for a tiu electron coupled to an Hg mode 

as a function of coupling. At weak coupling (dot-dashed), the ground-state energy 
(solid) drops 2.5 times more rapidly than the adiabatic minimum lowering (dotted) 
in order to reach the correct strong-coupling semiclassical limit (dashed) , which takes 
into account the softening of two (for n = 1) harmonic modes for a total zero-point 
gain hu! (indicated by a vertical arrow). (Adapted from Ref. [78].) 



is 5/2 times larger in perturbation theory than the simple adiabatic lowering [79]. 
This point is illustrated in Fig. 6.3 for n = 1. For a simple nondegenerate non-JT 
coupling, such as that to the Ag mode, the two quantities would instead coincide. 
The physical origin of this enhancement is a substantial reduction in zero-point 
energy. The ground-state energy drops faster at small g because the system trans- 
forms rapidly from 5 harmonic oscillators (zero-point energy = ^Huj) at g = to a 
Mexican- hat potential well, which is more "square- well" -like than the original hax- 
monic potential, and, in particular, carries several soft ( "pseudorotational" ) modes 
along the trough. Perturbative expressions for the ground-state energy lowering at 
weak coupling for ii" in all possible spin states are collected in Table 6.1. The 
perturbative expressions are essentially exact for g < 0.3. For larger (intermediate) 
coupling, the actual JT energy lowering is located between the perturbative "an- 
tiadiabatic" estimate (which systematically overestimates the JT energy lowering) 
and the adiabatic energy lowering (which consistently underestimates the energy 
lowering). 

In the strong-coupling limit, E^~^ concides with the classical adiabatic energy 

lowering (proportional to g'^hco), with the addition of a negative zero-point softening 
contribution (of order hu)). As was the case at weak coupling, in multi-electron 
contexts, JT distortions leads to a larger energy gain when charge concentrates 
in as few orbital as possible, thus favoring laxger distortions. Due to this, low-spin 
configurations end generally lower in energy than high spin ones, giving rise to a kind 
of reversed (first) Hund's rule exchange behavior (here we ignore the actual Hund 
rule, which is however important, and will be introduced in Sec. 6.2.4). To describe 
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Table 6.1. Perturbative (weak coupling, to order [23]) and semiclassical (strong 

coupling, to order gr"^ [72]) ground-state energies for an Hg mode coupled coupled 
to n electrons in a t\u orbital, for the allowed values of the total spin S. Here, 
only electron-phonon coupling is considered, while electron-electron repulsion is ne- 
glected. 
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correctly the quantum dynamics in the JT trough, a detailed analysis of the coupled 
electron-phonon system is needed. A crucial ingredient in this problem is the value 
or TT of an electronic Berry phase [22, 24, 80, 81, 82] 7 in the semiclassic motion: it 
imposes the boundary condition on the soft-mode pseudorotational dynamics in the 
trough, thus selecting vibronic states of a given parity oxp(i7) [22, 75], not unlike 
similar vibronic contexts [83, 84, 85, 86]. This boundary condition determines the 
leading corrections to the adiabatic energetics, of order hui/g^ , reported in Table 6.1 
[72]. 

The model Hamiltonian (6.2) is readily extended to include all vibrational modes 
of real Ceo. Each mode adds one harmonic term of the type (6.3) and a coupling 
term of the type (6.4), characterized by an individual coupling amplitude gk and 
frequency Wfc. In the evaluation of the total electron-phonon coupling, one should 
also include two Ag vibrons which also couple linearly to the LUMO electrons, even 
though they do not split its degeneracy. This is a simple polaron problem, which is 
exactly solvable [87] and independent from the tiu ® 8Hg problem. The amount of 
ylg-related electron-phonon energy is simply —n^g'^Tiuj/2, independent of the total 
spin S. 

Coming to the eight Hg modes, generally speaking, even within the harmonic- 
phonon linear-coupling approximation, a realistic description of the dynamical JT 
state of a CgJ" ion is a significantly more complicated affair than the single-mode 
problem. Linear superposition of in the coupling of individual Hg modes to a same 
tiu orbital is valid only in the antiadiabatic perturbative regime, but does not seri- 
ously apply to the coupling range generally accepted for Ceo: for a general coupling 
strength, there is in fact no linear superposition. The JT splitting of the LUMO 
increases with the cooperative action of the single modes, so that even moderate 
coupling of many individual modes builds up a rather strong total coupling. 

For the detailed values of the electron-phonon couplings, several estimates are 
available in the literature [33, 35, 36, 70, 88, 89, 90, 91, 92, 93], some of which 
are reported in Table 6.2. The main observation here is that individual mode cou- 
plings are indeed intermediate to weak. However, a significant spread in the values 
of different estimates is due to the large sensitivity of the coupling to the individ- 
ual phonon eigenvector, which differ largely in different approximations. Moreover, 
the total coupling of all calculations is significantly weaker than all experimental 
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Table 6.2. Dimensionless couplings gu for the eight Hg modes to the LUMO of Ceo, 
according to several theoretical and experimental determinations. In the literature, 
the electron-phonon coupling are often gauged in terms of the weak-coupling super- 
conductivity parameter A: the relation with the Hamiltonian coupling strengths gu is 
[88]: g^. = ^Xk/N{eF)/fiLUk for the Hg modes (where N{eF) is the density of states per 
Ceo at the Fermi energy ep; one would get gl = 3Xk/N{€F)/?iu)k for the Ag modes). 
9cB ~ Efcfffc measures the total interaction strength, and w^tf = X^/t 3fc'"'fc/ffcff 
provides and average phonon frequency to estimate zero-pint effects. The last line 
reports the total A/iV(eF) from Hg modes, according to the different estimates. 
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estimates, such as those obtained by fitting the vibronic features of the photoe- 
mission spectrum of 0^, [93]. This is rather surprising in view of the fact that a 
recent parameter-free calculation [94, 95, 96] found instead good accord with ob- 
served vibronic structures in photoemission data of neutral Ceo [97, 98], based on 
HOMO couplings computed [70] with the same ab-initio techniques as those for the 
LUMO reported in Table 6.2. Even though some ideas are being pursued to under- 
stand the theory-experiment discrepancy for the LUMO couplings [99], this is still 
an unresolved issue. 

As pointed out by Bergomi and Jolicoeur [100], experiments on anions in matrix 
may yield relevant information to the vibronic effects. Near-infrared and optical 
spectra of C"(7 ions in solution and frozen Ar are available [45, 101]. A major tiu 
tig optical transition near 1 eV is present for all values of n. It is accompanied by 
additional vibronic shake-up structures, typically near 350, 750, 1400 and 1600 cm^^, 
which involve the vibronic couplings of both the tiu LUMO and tig LUMO-l-1. 
Also, a satisfactory vibronic assignement of a recently observed gas-phase spectrum 
[39] is not yet available. The experimental information seems insufheiont as yet for 
any relevant comparison with our calculations. Well-defined vibrational spectra are 
instead available for chemisorbed Cg,^ [^9] and for A^Ceo alkali fuUerides [54]. In 
this case, however, interaction of the electronic tiu level with surface states or with 
other tiu states of neighboring balls turns the level into a broad band, and our 
treatment as it stands is invalid. Rapid electron hopping from a molecule to another 
interferes substantially with the dynamical JT process, in a way which is not known 



12 



Nicola Manini and Erio Tosatti 



in detail [76, 102, 103]. The spectra of negatively charged Ceo adsorbates and solids, 
in any case, do not present any clear evidence of vibronic splittings, but rather of 
the gradual continuous shift most likely due to a gradual overall change of geometry, 
also suggested by DFT-LDA calculations [37]. 

The full quantum-mechanical problem of a threefold degenerate electronic state 
linearly coupled to eight fivefold degenerate harmonic oscillators is conceptually sim- 
ple but in practice, for arbitrary couplings, rather hard to solve exactly. On one side, 
it is straightforward to compute the Hamiltonian matrix elements on a truncated os- 
cillator basis. Also, this matrix is sparse, involving nonzero terms only between states 
whose numbers of vibrons differ by exactly one. Moreover, the inclusion of states 
with up to N vibrons, because of linear e-v coupling, implies decay as exp(— AT) of 
components with A-vibron states, for large enough A. For this reason, a truncated 
basis set including all states up to A™^'' vibrons gives a variational estimate of the 
lowest eigenvalues, with good convergence with increasing iV™^ [22, 25, 95]. The 
structure of the problem makes it especially suitable for a Lanczos algorithm. By 
this numerical technique, not only the ground-state energy but also several low-lying 
excitations can be computed [104]. 

To get an estimate of the accuracy of the adiabatic and antiadiabatic approxi- 
mations applied to CgQ, the couplings of Ref. [70] (see Table 6.2) can be plugged in 
the strong and weak-coupling formulas of Table 6.1 and obtain the following esti- 
mates for the ground-state energy: E°^\ = —38 meV, -Bad+zp i ~ —131 meV, and 
^airtlad 1 — nieV. These estimates should be compared to the exact ground-state 
energy lowering obtained by exact Lanczos diagonalization: £|~act i — rneV 
[105]. The better accuracy of the perturbative antiadiabatic estimate is due to the 
relatively small-to-intermediate value of the total X]f;5fc = O.SA (Fig. 6.3 indicates 
that the zero-point corrected adiabatic energy J?^d+zp i becomes fairly accurate only 
beyond > 1.5). For n = 2 S = 0, and even for n = 1 if the larger couplings 

from PES [93] are taken, the adiabatic estimate comes in much better agreement, 
especially after a zero-point softening correction is included [105, 106]. In effect, this 
adiabatic-plus-zero-point-corrcction approximation should be quite generally appli- 
cable to all strongly coupled JT problems. The relative difference of the classical 
adiabatic energy from the exact energy shows the relative importance of phonon 
quantum effects in Cg,^ ■ So, while for n — 1 a weaker coupling case - the exact 
energy lies closer to the antiadiabatic expression, the exact result for n = 3 and 
even more for n = 2 is much closer to the adiabatic (zero-point corrected) estimate, 
because in CgQ and CgQ JT coupling is effectively stronger (see Table 6.1). 

Even with the uncertainties on the exact value of the couplings, fuUerene ions 
are generally taken as intermediate- to strong-coupled JT systems, the relatively 
strong coupling realized as the the cooperative effect of several moderate individual 
couplings. This fact is more general: any JT system with a large enough number 
of vibron modes weakly-coupled to the same electronic level behaves as if strongly 
coupled, since the contributions to the splitting of the electronic multiplet add co- 
operatively. In fact, for intermediate coupling, the JT energy gains of the individual 
modes do not add algebraically, as they would in the limiting perturbative and adi- 
abatic regimes. To further understand this, we may wish to see the splitting as if 
effectively due to the coupling to a single mode. In this picture, the addition of more 
coupled modes effectively acts like an increase in g, thus a rightward shift along the 
solid line of Fig. 6.3. Since the curvature is upward, the sum of the individual energy 
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Table 6.3. Perturbative (weak coupling, to order g^) and semiclassical (strong cou- 
pling, to order g~^) clcctron-phonon pair energies for an Hg mode coupled coupled 
to an average n electrons in a tiu orbital, assuming for all n a low-spin ground state. 
The rmmerical pairing energies are computed for the Cg,j couijliugs of Ref. [70] 
(see Table 6.2) in the anti-adiabatic approximation. The last column includes the 
contribution — ^ g'ifiuik = —93 meV of the Ag modes [107]. The strong-coupling ex- 
pressions are rigorously valid only in the limit where the intra-shell electron-electron 
repulsion can be totally neglected. 



n 


antiadiab 


adiabatic 

[9 00] 




[meV] 
[Hg modes only] 
(antiadiab) 


[meV] 
[Ag and Hg modes] 


1 

2 
3 




-g' + i- 




-195 
390 
-195 


-288 
297 
-288 



gains for each mode coupled separately is always larger than what one obtains from 
diagonalizing with all the modes included together. 

In a solid-state context, electron hopping between molecules will alter this pic- 
ture as soon as it is large enough to disturb substantially the local JT physics. 
However, as we will argue later on, the relevant hopping is not so much the bare 
electron hopping, but rather the quasiparticle effective hopping, a quantity which 
becomes much smaller near a Mott transition. So long as that effective hopping is 
small enough, the JT energetics as described above is expected to hold even in the 
solid state. Electronic correlation in a compound where conduction occurs through 
electron hopping between Ceo ionic sites is mainly governed by the extra electron- 
electron repulsion when the local charge fluctuates away from its average filling of 
n electrons. The pair energy, or Hubbard ?7, defined by 

IA„ = En+l + En-i — 2En , (6-5) 

(where En are the fully relaxed ground-state energies of n electrons) measures pre- 
cisely the strength of this local correlation. As we shall sec in detail below (Sec. 6.2.4), 
tin is largely dominated by the positive Coulomb electron-electron repulsion. How- 
ever, also the JT electron-phonon energies E^^ of Table 6.1 show a strong nonlinear 
n-dependence, and thus necessarily also contributes to W„. 

In the standard many-body perturbation theory language, formally lAn~^ is the 
real part of the two-electron vertex function at zero frequency, including electron- 
phonon contributions in the weak-coupling antiadiabatic limit. The contribution of 
a nondcgenerate Ag mode hCi ^ " = —g^huj < 0, independent of n and irrespective 
of spin S. If coupling to Ag modes prevailed, the electron-phonon pair energy would 
be negative, and thus electrons would lower their energy if, rather than keeping an 
uniform occupancy n of all sites, they separated into (n — 1) and (n + 1) at different 
molecules, thus creating a bipolaron, or a charge density wave. Even though charge 
disproportionation has been claimed to play some role in KsCeo [108], in practice 
the large value of the on-site Coulomb repulsive U parte rules this possibility out in 
the fullerides. 

Assuming a low-spin ground state for all n (wc> shall discuss the validity of 
this assumption in Sec. 6.2.4), Table 6.1 provides all the ingredients necessary to 
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Fig. 6.4. Exact single-mode electron-phonon pair-binding energy (solid line), 

compared to weak-coupling perturbation theory for <; ^ 1 (dotted) and semiclassical 
theory for g ^ 1 (dashed). Un~^ is similar for n = 1,3, 5 electrons. The coupling 
strengths for the Ceo modes are all less than ~ 1. (Adapted from Ref. [78].) 

determine approximate expressions for the Hg modes contribution to U^^^ , both in 
the antiadiabatic limit and in the adiabatic limit (neglecting retardation completely) . 
These expressions are reported in Table 6.3. For odd n, is negative, thus as 

customary the electron-phonon coupling is attractive and favors pairing [109], thus 
opposing Coulomb repulsion. Interestingly, for even n = 2, 4, U^~^ is instead positive: 
here the electron-phonon coupling of the JT type re-enforces the Coulomb repulsion 
and contributes to suppress local fluctuations away from the average number n of 
electrons at each site. The sign change of the JT contribution to U„ is due to the 
nonparabolic dependence of the ground-state energy on n, as reported in Table 6.1. 
Figure 6.4 shows a comparison of the approximate expressions with the pair energy 
obtained by exact diagonalization for n = 3. As for the total energy, the exact pair 
energy is bracketed between the antiadiabatic estimate (overestimating) and the 
semiclassical adiabatic expression (underestimating). 

The couplings of Table 6.2 can be plugged into the expressions of Table 6.3, to 
obtain an estimate of the phonon contributions to the pair energy. In Table 6.3 we 
use the (probably underestimated) couplings of Ref. [70], and the antiadiabatic ap- 
proximation (which insteads overestimates U„), to obtain our best realistic estimate 
of the phonon contribution to the pair energies. These values are relatively large, 
and even though they come a long way from reversing the much larger Coulomb 
repulsion addressed in Sec. 6.2.4 below, they represent a serious correction with 
cannot be neglected. 
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6.2.4 Intramolecular Coulomb repulsion and Hund's rule exchange 

We come now to the Coulomb electron-electron intra-molecular energetics. As ap- 
propriate to an orbitally degenerate level, two distinct Coulomb enter: the "Hubbard 
U" term coupling to the total charge in the degenerate shell, and an intramolecular 
exchange term. The first term is strictly related to the pair energy of Eq. (6.5): 
basically measuring the energy cost required to change electron occupancy of n of 
the molecular site. In a solid conductor, this term will oppose all local fluctuations 
away from the average occupancy. Exchange terms, of smaller value, act instead to 
split multiplets at any given fixed occupancy. They enforce in particular the first two 
Hund rules, lowering in energy the otherwise degenerate multiplet states of highest 
spin and largest orbital momentum. 

The structure of the Coulomb Hamiltonian for a (tiu)" configuration is formally 
the same as that for an atomic p" in spherical symmetry, and as such entirely 
determined by two parameters only, usually chosen as the configuration-averaged U 
and a Hund-rule intra-molecular exchange J [72, 110, 111]. For a detailed symmetry 
analysis of the structure of the Coulomb Hamiltonian H^^'^, we refer to Rcf. [112]. 

The Coulomb Hamiltonian is a 2-body electronic operator which can be written 

as 



are expressed on the basis of the single-particle LUMO wave functions ipmai^), 
associated to the fermion operators cl^„ introduced above. 

These Coulomb integrals (6.7) can be evaluated directly for the simple kernel 
Wcr,CT'(r,r') ~ ge/(47reo|r — r'l), with phenomonological molecular orbitals [113]. This 
approach neglects the screening due to the other electrons within the same fullerene 
molecule. Due to that, it yields very large exchange values, close to early Hartree- 
Fock results [114], which generally overestimate exchange duo precisely to lack of 
screening. With that in mind, one can regard as substantially more realistic the 
evaluation of the interaction Hamiltonian (6.6) obtained in Rcf. [112] on the basis of 
DFT-LDA electronic structure calculations, which account for the full polarization 
response of the total charge density. Since on the other hand LDA overestimates 
screening, the outcome of this calculation represents in turn an underestimate of 
actual intramolecular exchange. 

Like in the atomic case, the symmetry of the Coulomb interaction plus the 
molecular symmetry of the problem allow us to express all of the Coulomb integrals 
in (6.7) as functions of a small set of physical parameters. Indeed, the Coulomb 
integrals are expressed as follows 




(6.6) 



where the Coulomb integrals 




(6.7) 




(6.8) 



in terms of a minimal set of independent parameters F^, defined as: 
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= PjE l'^'^ /dV'<(r)u(r,r')<(r') 



(6.9) 



with 



<(r) = J2 Vm{v) ^n(r) . (6.10) 



The two parameters governing are related to the fc = and k = 2 Slater- 

Condon integrals F'*^' for p electrons in spherical symmetry [115]. 

With the decompositiou (6.8) in hand, it is convenient to re-organize the inter- 
action Hamiltonian (6.6) in terms of number-conserving symmetry-adapted fermion 
operators: 

H'-^ = \Y.F^{Y^ w^^ w^'n - i (f^« - 5F««) h (6.11) 

A \ IJ. J 

where we defined the operators 

(6.12) 



w 



Rather than the F parameters, it is customary to use the parameters U = 
pAg i^_pHg J ^ pHg 12 (^^jjg uotatiou K is sometimes used for this quantity; 

it also equals F^^^ /3 of spherical symmetry [72, 115]). With this definition, observe 
that 

E"-" ""''(n) = TrU(H°-°) = U ^ , (6.13) 

where Trjn is the trace restricted to the n-electrons states. U governs therefore the 
parabolic dependence of the Coulomb energy averaged over all possible multiplet 
configurations on the total number of electrons n. It should be noted that U in- 
volves average multiplet energies, not ground-state energies like the pair binding 
energy of Eq. (6.5): the two quantites therefore would only coincide in a hypotheti- 
cal system where multiplets were unsplit by JT and exchange terms. With standard 
multiplet splitting instead, the pair-energy definition, based on the ground-state en- 
ergy, depends on n, which makes it rather inconvenient to characterize the Coulomb 
repulsion in the tiu shell with a single microscopical parameter. 

The J parameter controls the multiplet exchange splittings, i.e. it implements 
Hund's rules. One can label all multiplet states with an orbital "angular momentum" 
L (recall that the t\u orbitals behave effectively as p orbitals), as in Table 6.4. In 
fact, it is possible to express the Coulomb energy of all multiplets with the equation 



~(n,S,L) = U ^ ^ ' - J 



2S{S + 1) + ^L{L + 1) + ^n(n - 6) 



(6.14) 



as a function of n, S and L [112]. These energies are also reported in Table 6.4. 

The Coulomb parameters were computed by ab-initio density-functional cal- 
culations in Ref. [112]. The value of the molecular Hubbard U term of isolated 
Ceo was determined U = 3.07 eV, in the same range as previous calculations 
[110, 116, 117, 118]. This estimate can be compared to experiment. The pair en- 
ergy Ui = Eo + E2~ 2Ei = (Eo - El) - {Ei - E2) = Ea{C6o) - Ea{C^o)' where 
Ea indicates the electron affinity. These quantities are experimentally available for 
molecular Ceo [40] and Cgg [43, 119] 
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Table 6.4. Coulomb energies for the (<iu)" or multiplets expressed in terms 

of the two Coulomb parameters U and J. These multiplets arc correct under the 
assumption that olectron-plionon coupling can be totally neglected. The n = 4 and 
n = 5 multiplet have the same J terms as n = 2 and n = 1, with the center-of-mass 
(U) term given by Eq. (6.13): 6U and 10 U respectively. 

eg- t^-'''+i)[L] t^^'+^U E---{n,S,Lf^ 
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Table 6.5. The energy (in meV) of the lowest state of for each n and S, including 
the electron-electron and JT contributions from Ho + He-v + H''~^ (but excluding 
the [en + Un{n — l)/2] term). Adiabatic approximation: the phonons are treated 
in the adiabatic (strong-coupling) approximation, by full relaxation of the phonon 
modes to the optimal classical JT distortion for each n and S. Anti-adiabatic: the 
electron-phonon behaves effectively as a negative J""^ = — |X)fc5fc^fc [112]: the 
combined Hund and JT interaction yields levels split according to the pattern of 
Table 6.4, with J replaced by J -h J"'"" . (Prom Ref. [112].) 
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adiabatic 


anti-adiabatic 
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-100 
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+75 



Ea (Cso) - Ea (Cgo) = 2.7eV - 0.17eV w 2.5 eV. (6.15) 

This value is slightly reduced from pure Coulomb by electron-phonon coupling and 
orbital relaxation (see also Table 6.3). If we neglect these comparatively small con- 
tributions to the pair energy Wi w J7: a range U ~ 2.5-^3 eV provides a realistic 
estimate of the actual U parameter of gas-phase Ceo- 

For the intra-orbital exchange interaction, the calculation of Ref. [112] finds a 
value J = 32 meV. Calculations where screening is ignored or underestimated find 
much larger values, e.g. J = 110 meV (Hartree-Fock [114]) and J = 95 meV (cal- 
culation of the unscreened integrals in model orbitals [113]). An early calculation 
includiug screening by the strongly polarizable Ceo molecule found J = 25 meV 
[110]. Because DFT tends to overestimate screening, DFT values of J may be some- 
what underestimated. We believe that the actual Coulomb parameters of Ceo could 
lie somewhere in between the DFT couplings assumed here and the "bare" ones of 
Ref. [113, 120], but most likely closer to the DFT ones, due to the large polarizability 
of Ceo. An intermediate value of J ~ 50 meV is probably realistic. 
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It should be noted in this respect that suggestions advanced in earlier times 
that ovcrscrccning could reverse altogether the sign of this frozon-moloculo, purely 
electronic exchange J [109, 121, 122] appear incorrect. This is confirmed by recent 
purely electronic quantum Monte Carlo calculations [123] finding a spin gap between 
triplet (ground) and singlet (excited) states consistent with J ~ 54 meV. What is 
emerging as the correct picture is that only after inclusion of nuclear motion, and 
specifically of JT effects that arc strongly anti-cxchango, the effective J may change 
sign, from positive to negative. As the above results in fact show, the JT energetics is 
opposed to Hund's rule, and favor low-spin states against intra^molecular exchange, 
as evident from Table 6.1. The effects of the two interactions, electronic exchange 
and JT, tend to cancel. In fact the antiadiabatic expressions of the JT energies 
(Tabic 6.1 and Rof. [23]) arc consistent with the structure given by Eq. (6.14), if 
we replace J with a negative local exchange J""^ = — | glhojk- Therefore when 
the antiadiabatic approximation holds - for example when the JT coupling is weak 
and the effective electron hopping is much smaller than important Hg vibration 
frequencies - one can account simultaneously for both JT and Coulomb exchange 
by replacing J with an effective Jeff = J + J"~^ . The resulting energetics for Cgg" 
ions is summarized in Table 6.5. Likewise, an additional term J''~^ ^n{6 — n) must 
be introduced to account for the shift of the center of mass and contributes U''~^ = 
2 je-v _ 1 glhcok to the multiplct center-mass Coulomb U. However, as (7 3> 
J, the cancellation of U is much less effective than for the exchange term, and the 
phonon's contribution leaves a strongly repulsive total U. We should stress again 
that for general couplings and general electron hoppings the electron-vibration and 
electron-electron contributions do not add so simply, and must be treated in full 
[124, 125]. Nevertheless those approximate results are quite useful in providing a 
basic insight in this competition. 

The balance between Hund's rule exchange and JT is much less definite for 
Cqq, where, according to calculations [105], overall Coulomb exchange marginally 
prevails, leading to regular Hund's rule high-spin ground states. In that case, for 
example C^q should have a S = 1 ground state as an isolated ion, and might or 
might not retain high spin in solution [59] or in solid-state compounds [61, 62, 63, 64], 
depending on the environment. 

Experiments do confirm a substantial screening of J. Consider first Cg^ , present 
in insulating A4C60 (A=K, Rb). NMR [126, 127, 128] has identified an important 
singlot-triplct excitation around 100 -i- 140 meV. Taking this spin gap as represen- 
tative of the '^S— transition between the singlet molecular ground state and the 
lowest triplet state (separated by 5 Jeff according to Table 6.4), NMR data are consis- 
tent with a small effective value of Jofr — —20-; — 28 meV. Thus the clcctron-phonon 
coupling counteracts very effectively the Hund's rule exchange, and eventually pre- 
vails, even if marginally, in these electron-doped fuUerene compounds. As we will 
see below, this is the origin of the low-spin character of the Mott insulating phases 
(Mott JT insulators), and of (singlet, s-wave) superconductivity in fuUerides. EPR 
and NMR of discrete salts indicate that also Cg^ ions have a low-spin ground state 
[129, 130]. The case of CgQ is more intriguing. Here, it appears that the balance 
between Coulomb effects, matrix effects, and JT energies is more critical, and the 
final ground state can jump between and '^S, depending on marginal effects. 
Early reports suggested that when frozen is in a triplet 

state [44], but this evidence has been questioned [15]. In discrete salts, the singlet 
and triplet have been claimed to coincide closely [131, 132], whereas in solutions 
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(a) (b) 

Fig. 6.5. (a) Band structure of KeCeo- Energy is referred to the tiu band maximum. 
The solid-state density of states is compared to Gaussian-broadened single-molecule 
levels. (Prom Ref. [139].) (b) Band structure of KaCeo. Energy is referred to the 
Fermi level. (Prom Ref. [1].) 



the singlet appears to be stabilized by about 70 meV [15, 133, 134]. Evidence of a 
10 meV spin gap attributed to long-lived fluctuating CgQ in cubic CsCeo suggests 
almost complete compensation, with marginal prevalence of the singlet state [135] . 
On the other hand, ferromagnetism of TDAE-Ceo [136] has been interpreted [137] 
in terms of fluctuating S = 1 local CgQ sites. Structural evidence of a JT distortion 
is relatively rare. When available it is not especially clear in showing twice larger 
distortion in C^Q [138] relative to CgQ , which is what is expected for low spin. The 
only indirect evidence of a distortion in K4C60 is provided by the splitting of the 
highest Tiu vibrational mode demonstrated in Ref. [7]. 

In summary, experimental and theoretical evidence agrees that (i) in the energy 
comparison between different number occupancies n (pair energy W„) the repulsive 
Coulomb energies dominate; (ii) in the energy comparison between different multi- 
plet (spin) states at fixed occupancy n (exchange term J), the attractive anti-Hund 
JT energies tends to prevail over slightly weaker competing Coulomb exchange, but 
there are exceptions and high spin may prevail in particular environments. 



6.3 Strong correlation in solid alkali fuUerides: general 

After this long preliminary on the properties of molecular ions we now come to the 
solid state compounds, that constitute our main concern. The starting point for 
the discussion of their electronic behavior is as always the standard electron band 
structure. The band structure of solid alkali fullerides has been computed by many 
authors [1, 3, 139, 140, 141], generally based on DFT-LDA methods. As ilhistrated 
in Pigs. 6.5 and 6.6, the bandwidth W of the LUMO band is generally obtained in 



20 



Nicola Manini and Erio Tosatti 





0.4 




0.2 




0. 


0) 








■0.2 




-0.4 



(q) 



(b) m.K^c^, 




X W X' 



5 10 

states / eV/ spin/ Cqq 



0.4 
0.2 



"S 0. 



-0.2 
-0.4 



(c) 



X W X 




5 10 

states/ eV/ spin/ Cjq 



Fig. 6.6. LUMO band structure of NH3K3C60 (a), obtained by an DFT-LDA 

calculation, compared to that of KaCeo (c). The corresponding densities of states 
arc given in panels (b) and (d). In panel (b) the density of states of K3C(io (dashed 
line) is re-drawn for direct comparison with that of NH3K3C60 (solid line). (From 
Ref. [140].) 



the 500 -i- 600 mcV region. Differences for different lattice structures are of course 
present, as well as a strong dependence upon the relative rotational ordering of the 
Ceo balls [142]. In particular, smaller bandwidths are obtained in correspondence of 
larger lattice spacings, associated, e.g., to larger intercalated cations such as Cs"*" in 
the place of K+. The main observation is that the electron bandwidths in the solid 
are generally a fraction of the separation between neighboring molecular levels, the 
latter giving rise to well defined rigid bands. In particular, the threefold LUMO 
band remains well separated from the LUMO+1 and HOMO bands through the 
Brillouin zone. The metal-insulator transition in the Cg,7 -based compounds can 
then be interpreted in terms of the localization of the electrons within the partly 
filled tiu LUMO band. Self-energy effects are predicted to increase W by about 30% 
in the usually more accurate GW approximation over that calculated by LDA [142]. 

A direct experimental determination of the bands in solid Ceo compounds is 
complicated by characteristic strong electron-phonon satellites derived from the vi- 
bronic couplings of Sec. 6.2.3, which tend to broaden and influence all spectra. 
Angle-resolved photoemission, which could access direcly the band dispersion is 
made additionally hard by the orientational disorder and the small size of the Bril- 
louin zone. A recent experiment [143] accessed directly the dispersion in a K3C60 
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Fig. 6.7. DMFT zcro-tcmpcraturc phase diagram for the two-band Hubbard model 
as a function of filling and of the ratio U/W of the local Coulomb repulsion to the half 
width. In this calculation molecular Coulomb exchange J = 0, and the possibilities of 
magnetism and/or charge order are ignored. The T = Mott transitions at integer 
filling are located between the pairs of finite-temperature estimates indicated by 
squares and diamonds. (Adapted from Ref. [152].) 

monolayer, obtaining an estimate W w 250 meV, significantly smaller than the 
DFT-LDA bandwidth of the monolayer. Newer experiments on KsCeo multilayers 
find an even smaller W ~ 150 meV, about one third of the DFT bandwidth [144]. 
These small values do not deny the DFT estimates but most likely represent a quasi- 
particle effective dispersion, whose magnitude relative to the bare band dispersion 
is strongly affected, and nontrivially reduced by correlations. It should bo noted 
that, even if there is so far no accurate prediction for the full fc-dcpendent spectral 
function, this kind of observed electron dispersion is not to be identified with the 
simple quasiparticle bandwidth zW either. The effective observed bandwidth prob- 
ably corresponds to the intermediate energy scale discussed in recent DMFT studies 
of a simplified model, where it is designated as T+ [145, 146, 147]. Indeed, the same 
angle-resolved photoemission measurement carried out on KgCeo, where correlations 
play no role, find [144] a band dispersion in very good quantitative agreement with 
the DFT calculation of Fig. 6.5 {W ~ 0.6 eV) [139]. 

In the light of the discussion of Sec. 6.2.4, it is clear that the largest energy 
scale of for the LUMO-band electrons in the fuUerides is the Coulomb repulsion U. 
It remains the largest energy scale, even with the extra screening characteristic of 
the solid state, which appears to reduce the molecular U from about 2.5 ^ 3 eV to 
a smaller 0.9 1.6 eV [116, 148, 149, 150, 151]. The smallest value in this range 
is probably closer to the correct bulk estimate in the metallic doped materials, 
while the upper value better characterizes insulating states and the molecules at 
the surface, where screening may be less effective due to reduced coordination [16, 
116]. The ratio of the local Coulomb repulsion to the bandwidth can therefore be 
estimated in the range 1.5 < U/W < 3 for typical LUMO-band electrons in standard 
fuUerides. 

A local Coulomb repulsion larger than the bandwidth is the signature- of strong 
electron correlation in a solid. Situations of this kind are conventionally described 
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Fig. 6.8. The critical Uc/W as a function of orbital degeneracy d at half filling 
n = d. For d = 1 band U^-^'"^'^ « 1.46 W. In these DMFT calculations the local 
Coulomb exchange is neglected (J = 0). The possibility of magnetic and/or charge 
and/or superconductive order is also ignored. (Adapted from Ref. [154].) 



in the language of Hubbard models. The single-band half-filled 3D Hubbard model 
is believed to undergo a metal-insulator transition of the Mott type, as soon as the 
ratio U/W exceeds some critical Uc/W of order unity [153]. The conventional small- 
U metal transforms to a largo-?7 Mott-insulator where electrons remain essentially 
frozen, one at each site, so that the occurrence of zero and double occupancies is 
strongly even if not totally suppressed. The residual effect of intersite hopping is 
to induce a generally antiferromagnetic correlation between the spins of electron 
at neighboring sites. Away from half filling the model has a metallic (and possibly 
superconducting) ground state even for large U/W. 

The single-band Hubbaxd model, however, can hardly be applied to the fuUerides, 
where there are d = 3 relevant bands, derived from the three degenerate molecular 
LUMO orbitals. The rough physics of the multi-band Hubbard model is summarized 
in the phase diagram of Fig. 6.7, for the simpler case of only d = 2 bands. At 
integer-filling, lines of strongly correlated Mott insulators extending above Uc/W 
are surrounded by the metallic phase encompassing all noninteger fillings. The same 
qualitative picture should hold for d = 3, where metal-insulator transitions should 
occur for stoichiometric phases, in particular n = 3 (as in KsCso) and n = 4 (as in 

K4C6o). 

Recently, the metal-insulator transition has been attacked theoretically within 
the dynamical mean field theory (DMFT) [155]. This method, which becomes exact 
in the limit of large space dimensions, has provided quantitative estimates of the 
Uc/W under various conditions. The lines of Fig. 6.7 arc constructed based on the 
DMFT [152], assuming a "pure" Mott transition, with no complications induced by 
possible spin or/and charge order in the insulating phase. The "paramagnetic" Mott- 
Hubbard transition occurs at different critical U for different orbital degeneracy d 
(number of bands) and fillings n. Earlier work showed that the Mott transition in 
orbitally d-dcgcncratc lattice models takes place at larger values Uc/W for larger 
degeneracy d, roughly proportionally to y/d [156, 154, 157], as illustrated in Fig. 6.8. 
Specifically, the Mott transition in the half-filled one-band model (d = 1, n = 1) 
has been calculated at Uc'^'"^'^/W « 1.3 in a previous DMFT study based on the 
the quantum Monte Carlo method at the relatively low temperature T = W/32 
[158]. A more accurate T — estimate of Uc~'^^^'^ /W « 1.46 is provided by iterated 
perturbation theory [155]. The transition was found at U^-'°'''"^/W w 1.5 for d = 2, 
n=l [152], and Uc~'^'''"^/W w 1.8 for d = 2, n=2 [152, 154], with the same method. 
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and at the same temperature. T = Lanczos diagonalization DMFT calculations 

push these values slightly up, with the n=l value close to t/^ "'"'"'^/VK ^ 1.8 [140]. 
For d = 3, n=3, the transition was located close to U^''"'''^'^ /W ^ 2.3 [154], ahhough 
at a slightly higher temperature. The finite temperature appears to affect somewhat 
the numerical values: the zero-temperature DMFT calculations based on Lanczos 
diagonalization provide slightly larger values [159], but confirm that [/^-band ^ 

These values of Uc/W are significantly reduced when more realistic DMFT cal- 
culations allow for intrarsite couplings and the ensuing magnetically/charge ordered 
phases, and when any kind of local exchange term is included to break the local mul- 
tiplct degeneracy. Therefore one cannot compare directly the values of Uc/W from 
paramagnetic DMFT calculations of the multi-band Hubbard model with the metal- 
insulator transition in the actual compound, where it is affected by the detailed band 
structure, nesting, JT electron-vibration interaction and intra-site exchange. Most 
of these effects tend to lower Uc. 

In fact, the range 1.5 < U/W < 3 typical of the alkali fuUerides should put 
all these materials above the metal-insulator transition: all integer-filled compounds 
should then be insulators. In practice, the alkali fuUerides lie experimentally close 
to the Mott transition, on both sides. Different compounds exibit a variety of be- 
havior, including unconventional metals like cubic CsCeo [135], superconductors of 
the A3C60 family (A= K,Rb,Cs) [13, 151] including also NH3Na2CsC6() [160], a fer- 
romagnet - TDAE-Ceo [50, 51, 136, 137, 161, 162], and Mott insulators. Among the 
latter we can place A4C60 [4], possibly Na2C6o [163], (NH3)6Li3C6o [12], and the 
full class of compounds NH3 K3C60, NH3 K2RbC6o, NH3 KRb2C6o, and NH3 Rb3C6o 
[6, 8, 9]. The fact that the n = 2, 4 compounds are insulators is not expecially sur- 
prising in view of the increase in repulsive pair energy due to olectron-phonon con- 
tribution (Sec. 6.2.3), and the probably smaller Ui^-^'"^'^{n = 2) < U't^'"^'^{n = 3). 
Altogether, this scenario indicates a reduction of U in solid fuUerides, and a lowering 
of Uc relative to its bare theoretical value of Figs. 6.7 and 6.8. 

Magnetic susceptibility data [12, 164] in particular place the metallic/superconducting 
CgQ compounds close to the metal-insulator transition. Indeed, in the Hubbard 
model, the strongly-correlated metallic state near Uc is expected to show an anoma- 
lously large uniform magnetic susceptibility x compared to the free-electron value 
/x|Af(ep) (where N{eF) is the density of states per Ceo at the Fermi energy ep) [155]. 
The magnetic susceptibility is observed to increase rapidly with the lattice spacing 
in the Cgg fuUerides [12, 164]. When moving closer to the metal-insulator transition, 
the observed increase cannot be interpreted purely in terms of band narrowing and 
a related increase in the density of states at the Fermi level N{e,-p)- The discrepancy 
between this increase due to bare band narrowing and observation is huge, and can 
only by explained by strong correlations. 

Further, if less direct, evidence of the strong electronic correlation in K3C60 
is provided by the electron spectral function at Ef which can be assimilated to 
a narrow, dispersionless Kondo-like resonance by photoemission data [165]. On the 
other hand, photoemission is on one hand strongly affected by vibrational efl^ects and 
also possibly more sensitive to surface molecules [166], where correlation is stronger 
due to a less screened Hubbaxd U, thus caution is required before photoemission 
data should be taken as representative of bulk KsCeo as magnetic susceptibility 
data can. 
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Fig. 6.9. Schematic electronic phase diagram of Cg(7 compounds, showing the ap- 
proximate location of the metal (superconductor)-insulator phase boundary. The 
open symbols are literature values of Tc for a variety of superconducting ful- 
lerides, while the solid symbols mark Tjv (ambient pressure) and Tc (> 1 GPa) 
of NDsKsCeo- This expanded trivalont fuUcride is only metallic and superconduct- 
ing under pressure, and reverts to an antiferromagnetic insulator at zero pressure, 
in strong analogy with the cupratcs (exchanging pressure with doping). Moreover, 
the CgQ is in a S = 1/2, low-spin state, and not in S = 3/2 high-spin state. This 
indicates an effective inversion of the effective Hund's rule J, supporting the con- 
clusion that ND3 K3C60 and related compounds are Mott-JT insulators. (Adapted 
from Ref. [167].) 

In the ammoniatcd compounds, insertion of the electronically inert NH3 molec- 
ules expands the Ceo lattice, turning the cubic, metallic, and superconducting state 
of KaCao into an orthorhombic narrow-gap S = 1/2 antiferromagnetic insulator 
[6, 8, 168]. Upon application of pressure to NH3K3C60, the Mott insulating state 
reverts to a fully metallic and superconducting - while still orthorhombic - phase 
[169, 170]. A similar pressure-induced reversion from insulator to metal has been 
observed in tetragonal Rb4C6o [126]: these experiments rule out the possibility that 
metallic states of the fullerides should entirely be attributed to nonstochiometry. 
These transitions in the alkali fullerides are rare examples of experimentally ac- 
cessible Mott transitions. We will take that of NH3 KsCeo as a paradigm where 
several of the concepts illustrated above emerge most clearly. It was suggested [167] 
that the increase in volume per Ceo molecule relative to K3Ceo, with its probable 
slight decrease of electronic effective bandwidth W, should drive the Mott transition 
(Fig. 6.9). 

While the above is certainly relevant, in addition to the simple change in volume, 
thus of W, a role of the crystal-structure anisotropy in this pressure-driven transition 
from the insulator to the metal and superconductor has been proposed [8]. Crystal 
anisotropy affects hopping between neighboring fullerene balls, making it different 
for pairs of molecules at different distances, along different crystal directions, thus 
affecting the band structure. The effect of anisotropy is most clearly seen in the the 
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band splitting at F, illustrated in Fig. 6.6. A strong enough "crystal-field" splitting 

of the threefold degenerate tiu molecular orbital of Ceo caused by orthorhombic 
lattice anisotropy could remove one or several bands away from the Fermi level, 
effectively reducing the orbital degeneracy d, and shifting the Mott transition to 
a smaller critical value Uc/W. Exploring this concept, the basic question is how 
large a splitting is required to promote the effective reduction of degeneracy. In a 
non-interacting system, that reduction would clearly require a splitting magnitude 
at least similar to the electron bandwidth W. In a strongly- interacting system, it 
is important to understand whether the effective degeneracy lifting, and the asso- 
ciated substantial displacement of the metal-insulator transition, will again require 
anisotropic splittings as large as the bandwidth, or else if some smaller energy scale 
will emerge in its place. 

In the final sections in this chapter we will adddress this issue within DMFT, 
by studying the effects of a band splitting on the Mott transition of an orbitally 
degenerate, strongly correlated metal. The theory for the distorted ammoniated 
compounds underscores the central role of a new low-energy scale characteristic of 
strongly correlated itinerant electrons close to the Mott transition: the width zW of 
the Kondo peak. We will axgue about the importance of this concept in the context 
of both CgQ and CgQ fuUerides, and of strongly correlated materials in general. 



6.4 Strong correlations in fullerides: theoretical models 

Models of alkali fullerides should be able to describe metals, Mott insulators, and su- 
perconductors. The largest energy parameter being the intra-molccular repulsion U, 
with the electron hopping only second, ab initio approaches are not yet in adequate 
shape to meet this challenge at the present date. 

Early work attempted to interpret the superconductivity of the fullerides in 
terms of Migdal-Ehashberg theory [22, 23, 88, 89, 171, 172, 173, 174, 175, 90]. The 
necessity to introduce a comparably large value of the pseudopotential fx* w 0.4 
already suggests that electronic correlations are exceedingly important. The large 
zero-point vibrational energy ^hui associated to the vibrationally-derived optical 
modes of these solids also suggests a breakdown of Migdal's approximation, and 
that vertex corrections and other nonadiabatic effects could play a relevant role in 
the fullerides [176, 177, 178, 179, 180, 181). However, this kind of theory starts from 
free electrons in a regular band metals, hard to accept for the fullerides whose large 
value of U puts them close to (if not beyond) the Mott-Hubbard transition. As an- 
ticipated above, for a satisfactory description of the large local Coulomb correlation 
and of the band physics on the same footing, a multi-band Hubbard model, possibly 
including the JT phonons, and described within the DMFT [155] or some of its 
several evolutions [182, 183, 184] is certainly a more natural choice. 

The simplest d-band Hubbard model is written as 

H ^ Hbop + H"-" . (6.16) 

Here we assume purely diagonal hoppings between orbitals at nearest neighbor sites 
i and j: 

d 

Hhop = {4ma^jrn.a + H.C. ^ + ^em flim, (6.17) 

,a m=l i m 
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with d bands of the same width W, the anisotropic symmetry lowering entirely 
embodied in a local diagonal term splitting the on-site orbital energy em- Of course, 
these simplifying assumption are not strictly applicable to the fuUerides, where a 
more complicated 3x3 matrix of overlaps for the tiu LUMO orbitals at neighboring 
sites (possibly including merohedral disorder) should be (and has been) employed 
[16, 185]. The essence of the Mott transition in the fuUerides should most likely 
not be affected qualitatively by the details of the lioppings, but the quantitative 
previsions of a model ignoring the correct tight-binding hopping matrix, including 
effects of orientational order or merohedral disorder, should not be taken too literally. 

The Coulomb part if""'* = He-ei, with the local term defined in Eq. (6.11) 
for d = 3. Another equivalent formulation for the local term goes as follows: 

+ {u-j) J2 

(6.18) 

m m<m' a 

+ U nimtfiim'i + J ^ '^Im't'^iml^im' iCimt , 

with the advantage that it applies equally well for any degeneracy d = 1, 2, 3 [154]. 

The on-site electron- vibration couplings Eq. (6.4) and the ensuing JT effect, 
are quite important for many aspects, including superconductivity and resistivity 
[16, 23, 151, 186], as discussed in Sec. 6.4.4 below. While we included it in the single- 
ion model of Eq. (6.16), the initial simplification of neglecting the JT coupling in 
the solid is useful for a simple approach to the basic physics of the Mott transition. 

6.4.1 The bare d = 2 model with anisotropy 

For an initial study of the anisotropy effect on the Mott transition, one can start 
with the simplest orbitally degenerate Hubbard model, namely d = 2 bands, zero 
JT coupling, zero Hund's rule exchange. For the anisotropy term we can assume, 
without loss of generality, £2 = — ei ~ A/2. We choose to study filling n = 1 
electron/site: this choice is motivated by the observation that for d = 2, half filling 
n = d = 2 would yield a trivial phase diagram, not comparable to the realistic case 
n = d = 3. We study this model at zero temperature, where the Mott transition 
appears most clearly, as a function of A/W and U/W. 

In the [7 = limit, the splitting A > simply shifts band 2 upwards with 
respect to band 1, promoting electron transfer from the upper to the lower band. 
Above a critical value A = Ac, the upper band will be emptied. For example, with 
two symmetric bands, above Ac/W = 0.5 the upper band is emptied and the lower 
band remains half filled. At ?7 = 0, the transition between the "two-band metal" and 
"one-band metal" is continuous. Because the topology of the Fermi surface changes, 
this transition is accompanied by a weak zero-temperature singularity of the total 
energy first described by Lifshitz [187]. When the electron-electron interaction U 
is turned on, one expects the emptying of the upper band to take place at smaller 
values of Ac, owing to the effective band narrowing. Perturbatively in U one can 
show that at weak coupling 

Ac{U) = Ac{0) - -r U + 0{U^) , (6.19) 

where the value of the coefficient 7 > depends on details of the bands. In addition, 
electron-electron interactions might modify the nature of the metal-metal transition 
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singularity relative to the noninteracting case [188], a point which we will not further 

address hero. 

The A = and A ^ oo limits reduce then, respectively, to the (quarter-filled) 
two-band and (half-filled) single-band Hubbard models, both possessing a metal- 
insulator transition as a function of U [152, 158, 159, 189]. We are not interested 
here in the weak-coupling antiferromagnetic instability of the ensuing large- half- 
filled band, associated with the possible presence or absence of nesting characteristic 
of a specific assumed intersite hopping scheme. Consistently with our neglect of all 
intrasite multiplet interactions, we also ignore for the moment the possibility of 
charge ordered and/or superconducting phases. In particular, we leave the local 
exchange term out (J — 0), with the proviso that this term will later be crucial in 
order to understand the ordered phases which enrich the phase diagram of the model 
[146, 186, 190, 191]. We assume a genuine Mott transition for the half-filled single- 
band model to occur at a finite U^'^'""^ > (for A ^ oo) and U^''^'""^ > U^-'^'""^ 
(for Z\ = 0). 

The limit of strong interaction, U ^ W , is insulating for any value of A. This 
limit can be studied by mapping the model (6.16) onto a spin and orbital exchange 
Hamiltonian which reads [137] 

-ffexch = JtY, (^i ■ +Ti-Tj+ASi-Sj Ti-Tj)-AY,Ti , (6.20) 

where the true spin operator is Sj = i Ximi/i/' '^]mi^ '^'^v' '^jmv' ^■^d the pseudospin- 
1/2 vector operators Tj = \Y,mm.'u'^]mv'^mm' Cj^'u represent the orbital degrees 
of freedom, a being the vector of Pauli matrices, and Jt = 2t^ /U . For Z\ = this 
model has been studied both in one [192] and two dimensions [193], with suggestions 
that interesting spin-liquid physics could be realized. For our purposes, it suffices 
to note that this model has no ferro-orbital instability, and has therefore a finite 
5 = orbital susceptibility. As a consequence it takes a nonzero value of A to fully 
orbitally polarize the ground state. Due to the absence of cross-band terms in the 
kinetic energy, complete orbital polarization occurs at a finite Ac (x Jt = 2t^ /U. For 
A > Ac the ground state is thus represented by a one-band Mott insulator plus a 
totally empty split-off band. 

For the specific case we arc interested in (n — 1 electron in a two-fold degenerate 
band) no orbital ordering is present at weak coupling. On the other hand, the pos- 
sibility of antiferro-orbital and/or spin ordering within the Mott insulating phase, 
and of spin-orbital density waves in the intermediate U regime, depends crucially on 
the details of the various hopping matrix elements and of lattice coordination. These 
are left out of the infinite-dimensional lattice assumed in the calculations to be de- 
scribed in the following section. Therefore, within the standard scheme of single-site 
DMFT calculations, we can study the phase diagram of this model restricted to spin 
and orbital paramagnetic states. 

Figure 6.10 is a sketch of the zero-temperature phase diagram of the d = 2 model 
in Eq. (6.16) as a function of {U,A) for n = 1 electron per site. The transitions at 
Zi = and A = oo are located at the values of U/W ~ 1.8 and 1.5 as discussed 
in Sec. 6.3. The descending AB line in Fig. 6.10, representing Uc as a function 
of A separating metals from Mott insulators, indicating that [Z^"*"""* > Uc'^'"^'^. 
The DM' line separates the fully orbitally polarized Mott insulator (on the right) 
from the two-band insulator, roughly as Uc ~ A~^ for small A. Similarly, the CM 
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Fig. 6.10. Qualitative zero-temperature phase diagram for the two-band Hubbard 
model at quarter filling (one electron per site) in the U-A plane, where A is the 
anisotropy splitting of the two orbitals. The various phases and lines are described 
in the text. The multicritical points M and M' are not necessarily distinct. (From 



line separates the fully orbitally polarized metal from the two-band metal: it starts 
from point C with a linear slope —7"^, according to Eq. (6.19). In the region of full 
orbital polarization the value of A is irrelevant, and this is the reason why the Mott- 
transition line MB is horizontal. As A increases, for U < f/c"*'""'* the upper-band 
emptying transition takes the two-band metal across the CM line over to a one-band 
metal, while for U > U^-^'"^'^ it leads across the AM line to a Mott insulating state. 
The steeply dropping AM line is the main outcome of the calculation, showing that 
the two-band Mott transition is heavily affected already at small A <C W, and not 
at Z\ ~ W, as one might have expected. As the DMFT calculations below show, 
the effective emptying transition occurs when A increases to reach Ac{U) oa zW (x 

Further quantitative informations concerning this phase diagram (in particular 
in the region of intermediate 17 ~ and A < W) were obtained [140] by means 
of DMFT in the exact diagonalization flavor [155, 159] on the Bethe lattice. The 
results of the DMFT calculations are summarized in Fig. 6.11. The DMFT critical 
fji-band ^ 1 and U^-^'''"^ ~ 1.8W are in fair agreement with corresponding 

values obtained by other methods [152, 158, 159, 189]. The other points in the phase 
diagram are obtained by following the stability of the two-band metal for a given 
value of U and increasing A, marking the emptying transition to the one-band metal 
or to the insulator. 

A deficiency of this single-site DMFT calculation - which is restricted as usual 
to paramagnetic states only - is the absence of a two-band insulating state for A > 
(i.e. the M' point coincides with A). In fact, the suppression of the antiferro orbital 
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Fig. 6.11. DMFT zero-temperature phase diagram for the two-band Hubbard 
model at quarter filling (one electron per site) , obtained by the exact diagonalization 
method in the paramagnetic sector. (From Ref. [140].) 



fluctuations embodied in the exchange model (6.20) produces a fictitious infinite 
uniform orbital susceptibility which leads to full orbital polarization as soon as A 
is turned on [155]. Despite this limitation, the results of the DMFT calculations are 
suggestive, revealing the announced sharp reduction of the metal-insulator Uc{A) for 
small but finite A. Indeed, the DMFT yields a Ac{U) which is roughly proportional 
to the quasiparticle residue z{U) of the undistorted {A = 0) correlated metal at 

(X l3z{U) , (6.21) 

with a proportionality constant f3 ~ 0.3. Since z{U) vanishes as (7 -> (Jc"''"-"'^, most 
likely linearly in {U^-^'''"^ - U) [155], even a small A <^W is sufficient to cause a 
metal-insulator transition in the strongly correlated metal. For example, following 
the bold arrow &t U = 1.5 W in Fig. 6.11, a A value as small as 0.08 W is sufficient 
to cross the transition line from the metal to insulator. 

Illustrating further the transition. Fig. 6.12 shows the behavior of the spectral 
density 

Gm{uj) being the one-particle Green's function of band m, on both sides of the 
metal-insulator transition. The asymmetry of the upper-band spectral density A2{uj) 
(solid lines) is very pronounced, as this band is nearly {A/W = 0.07) or completely 
{A/W = 0.08) empty. As soon as the Kondo-like peaks of the two bands differ 
enough in energy to induce the emptying of band 2, the lower-band spectral density 
Ai (lu) takes the symmetric shape, characteristic of the half-filled one-band Hubbard 
model. Here the quasiparticle peak disappears completely, as this value of (/ > 
^i-band p^^g ^j^g jjubbard model of band 1 well inside the insulating regime. 
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Fig. 6.12. Spectral density = — 7r~^Im G'„i('^) at \J jW = 1.5 across the 

Mott-Hubbard transition for increasing anisotropy splitting Z\. a; is referred to the 
Fermi energy. Solid lines refer to the minority orbital m = 2. The multi-peak struc- 
tures of the high-energy side bands are artifacts of the finite discretization of the 
conduction band. (From Ref. [140].) 

6.4.2 Noncubic band splitting and correlations in NH3K3C60 

The above d = 2 model calculations show that a small splitting A cc zW ai the or- 
bitally degenerate band can drive the metal-insulator transition. We wish to explore 
the implications that this result - if assumed to be more general than the simple 
model where it was derived - can have on the metal-insulator transition which takes 
place between isoelectronic K3C60 and NH3K3C60 (the former cubic and the lat- 
ter orthorhombic) and on the insulator-metal transition of NH3K3C60 itself under 
pressure. 

The effect of ammoniation is apparently twofold. The first effect is a large volume 
expansion, implying some band narrowing as well as an increase of U due to reduced 
screening. The second effect is a breaking of cubic symmetry in the anisotropic 
lattice structure. The actual strength of anisotropy in the ammoniated fuUeride 
may be estimated quantitatively from the DFT-LDA bands of Fig. 6.6. They were 
obtained in a simplified geometry with a single Ceo unit cell with lattice constant 
a = 14.2 A for fee K3C60, and a = 14.89 A for NH3 K3C60 the latter with a centered 
tetragonal unit cell with c/a = 0.91, neglecting the exceedingly small difference 
between a and h. Merohedral disorder and the rich antiferro-rotational structure 
recently discovered in actual NH3K3C60 [194] are also ignored. The bands indicate 
that the insertion of the NH3 molecules modifies only weakly the essentially pure 
K3C60 LUMO conduction band, as expected. In particular, as apparent from Fig. 6.6, 
the DFT-LDA bandwidths of NH3 K3C60 and K3C60 are very similar, both W ~ 
0.6 eV. The main difference in the two band structures is a splitting of the threefold- 
degenerate t\u band of K3C60 at the F point of NH3 K3C60. The ^T-point splitting. 
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Fig. 6.13. Simplified geometry of NH3 K3C60 (top) and K3C60 (bottom) used in the 
DFT-LDA calculation. Octahedral K atoms are indicated in black, while tetrahedral 
ones are in gray. All visible atoms of the central unit cell are shown together with 
some atoms in the neighbor unit cells. In NH3 K3C60 the c-axis is in the z direction. 
(From Ref. [140].) 

roughly corresponding to a dimensionless ratio A/W ~ 0.25, provides a measure 
of the strength of the non-cubic crystalline environment seen by the tiu orbital on 
each fuUerene molecule in the orthorhombic structure of NH3K3C60 . Its magnitude is 
roughly a quarter (0.2 — 0.3) of the total bandwidth; this represents the main result 
of this DFT calculation. Now if this compound were an uncorrelated metal, this 
splitting would have no real consequences, besides a change of shape of the Fermi 
surface. The consequences can be much more important due to strong correlations, 
as anticipated above. 

In fact, as discussed above, the fuUerides lie on the metallic side of the Mott 
transition in the half-filled d = 3 band cubic system. Upon ammoniation of K3C60 
into orthorhombic NH3 KaCeo, the volume expansion will first of all increase U/W, 
while the anisotropy will give rise to a nonzero A, corresponding to a displacement 
as indicated by a rightward arrow in Fig. 6.11. Inspection of Fig. 6.11 shows that, if 
U is close enough to Uc, even a band splitting A substantially smaller than what we 
have estimated for NH3 K3C60 could suffice to drive that metal insulator transition, 
even without any appreciable change in U/W. Conductivity measurements [195] 
on the class of compounds NII3 Ks-ajRbajCeo supports the possibility that the or- 
thorhombic distortion could be an important ingredient driving the Mott transition 
in these systems. An interesting - even if perhaps not practically straightforward 
- test of this overall picture could be obtained by applying uniaxial stress to cubic 
superconducting fuUerides of the A3 Ceo family. Contrary to the standard tendency 
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Fig. 6.14. A schematic (l],^) phase diagram for the d = 3-bands family of com- 
pounds KaCeo and NHaKaCeo- The dashed and dot-dashed arrows indicate two 
plausible paths of pressure-induced metallization of NHsKsCeo- Cubic insulating 
(NH3)6Li3C6o [12] is also indicated. (Adapted from Ref. [140].) 

of hydrostatic pressure toward metallization and lower superconducting Tc, the or- 
bital splitting associated with the appropriate uniaxial strain could drive some of 
these compounds towards stronger correlation, thus possibly first toward higher Tc, 
and eventually Mott insulating. A similar suggestion was put forward by Koch [196] , 
although on rather different grounds. 

We must stress here that the scenario sketched above does not at all diminish the 
role, in the metal-insulator transition, of the accompanying expansion of the lattice, 
and of the dependence of U /W upon volume. That role is experimentally proven by 
the observation of a pressure-driven insulator-metal transition, where NH3 K3C60 
is transformed into a metal (and a superconductor), despite the permanence of the 
orthorhombic structure [170]. 

6.4.3 Mott transition in other fullerides and other molecular 
conductors 

The vicinity to the metal-insulator transition is not exclusive of NH3 A3 Ceo- Other 
systems of the same family are (NH3)6Li3C6o [12], (NH3)a; NaRb2C6o {x — 1.6) 
[197], (NH3)a;NaK2C6o (0.5 < a; < 1) [198, 199, 200], and the noncubic expanded 
alkali fulleride CS3C60 [201]. The pressure-driven insulator-superconductor transition 
observed in the latter compound further underscores the role of volume expansion 
in favoring the Mott insulator. In fact, (NH3)6Li3C6o is an example of a compound 
of the n = 3 family lying beyond the Mott transition despite its cubic structure [12], 
as illustrated in Fig. 6.14. 

The n = 4 alkali fullerides also sit in the vicinity of the metal-insulator transi- 
tion, usually on the insulating side. Indeed, the pressure-induced insulator to metal 
transition observed in tetragonal Rb4C6o [126] is reminiscent of that observed in the 
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Fig. 6.15. Figure 1. Electronic phase diagram of organic conductors. The verti- 
cal axis indicates electron correlation U/W. The organic antiferromagnetic Mott 

insulating state (k-ET2Cu[N(CN)2]C1) transforms into a superconducting one [n- 
ET2Cu(NCS)2], with a reduction in U /W. The horizontal axis shows the band-filling 
control. (From Ref. [202].) 

Cg(7 compounds NH3 K3C60 and CsaCeo, but for the lack of superconductivity. We 
shall come back to these Cgg compounds in the final discussion. 

In the wider class of molecular conductors, the Mott transition is not exclusive of 
the fuUerides; other organic compounds show a similar phase diagram. Especially the 
k;-ET2 [ET = bis(ethylenedithio)tetrathiafulvalene] compounds show a similar phe- 
nomenology to trivalcnt fullcridcs [202] , with a clear Mott transition at stoichiomctry 
and a superconducting phase all around, Fig. 6.15. More recently, insulating metal- 
phthalocyanine organic films were turned conducting and metallic by alkali doping 
[203, 204]. A similar scheme to fuUerides, and to the metallic organics, has been 
proposed [205] as a possible scenario for the electron-doped metal-phthalocyanine 
conductors. The verification of the existence of stoichiometric phases, of Mott insu- 
lators at integer fillings, and possibly of superconducting phcises in between remains 
a challenge for future experiments and theory in these materials. 

Next, we should mention compounds such as (AsF6)2 Ceo [61, 62, 63, 64], where 
fullerene is stripped of two of its h„ valence electrons. The physics of positive Ceo 
ions has been described in Ref. [105], and is somewhat parallel to that of negative 
ions, though with a richer parameter space which implies several complications. The 
Hubbard U in the hu HOMO level is slightly larger than in tiu LUMO. Cancellation 
between Hund's rule exchange and JT effect, discussed in sections 6.2.4 and 6.4.4, is 
present here too. It was calculated to be nearly exact in the doubly positive ion, with 
a slight prevalence of exchange and a marginally stable high-spin ground state [105]. 
The doubly doped-solid state compounds (AsFe)2Ceo and (SbFe)2Ceo have been 
reported to be small gap insulators [61, 62, 63, 64]. There appears to be no simple 
way that they could be band insulators, and the likeliest explanation is that they arc 
Mott insulators. As for magnetism, some indirect experimental evidence has been 
suggested, consistent with high spin in the di-cation. However the balance between 
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intramolecular exchange and JT is very marginal, which leaves the possibility for 
a low-spin holc-dopcd Mott insulator wide open. There is oven some preliminary 
indication that this might in fact bo realized in (AsF6)2Cb() [206]. It would be very 
interesting in this respect to begin an experimental study of the insulator-metal 
transition in these compounds, to be obtained for example by hydrostatic pressure. 
Especially if the transition was relatively continuos and not too strongly first-order, 
the ensuing metallic phase close to the singlet Mott insulator should exhibit a pseu- 
dogap with strong pair correlations, and possibly superconductivity, similar to that 
described in Ref. [191]. 

Finally, we note that the pressure- or doping-induced transformation of Mott 
insulators into metals and strongly correlated superconductors which we discuss 
here for fuUercnc compounds, is a scenario which shares many important elements 
with the high- Tc cuprates. Again, this is an issue to which we will return later. 

6.4.4 The crucial role of exchange and of JT on-site interactions 

While the above bare Hubbard model defines the main backbone of the phase di- 
agram, the precise details depends strongly on a variety of secondary, on-site cou- 
plings, generally present as permitted by the orbital degeneracy of the site. They 
act as a rule such as to optimize the energetics of the isolated site. As such, they 
generally favor the Mott insulating state, where electrons effectively localize on site. 
In the metallic state, conversely, kinetic energy washes out much more the effect of 
these on-site couplings. Electrons in the insulating state where on-site occupancy is 
close to integer and poorly fluctuating, can take full advantage of these secondary 
couplings, whereas travelling electrons in the metallic state cannot. As a result the 
overall effect of secondary on-site couplings such as J and Ejt is to reduce the crit- 
ical Uc- On the contrary, we can expect that strongly nonlocal/intersite couplings 
could favor the metallic state, and raise Uc . In the fuUerides, intramolecular couplings 
dominate, while screening may be assumed to suppress the long-range interactions. 
In particular, while we certainly cannot rule out some role for an inter- molecular 
Coulomb interaction V , we will completely neglect it here. 

As has been repeated over and over, the two important intra-site interactions 
are the JT coupling, Eq. (6.4), and the Hund's rule exchange part of Eq. (6.11). In 
the DMFT calculation of Sec. 6.4.1 however, these terms were initially left out. In 
fact, these on-site interactions, even if small in comparison with U and W , play a 
crucial role on the Mott transition in a band-degenerate conductor. Their presence 
can actually change nearly everything: the nature of the Mott insulator, its spin 
and orbital contents, the nature of the metal phase, its tendency to superconduct 
or not, and of course the precise location Uc of the Mott transition itself. Probably 
the simplest manner to illustrate this physical situation is to describe the results of 
DMFT calculations with proper inclusion of on-site interactions. 

On a given molecular ion, Hund's rule exchange will favor maximum spin and 
maximum orbital momentum, while JT will favor just the opposite. Because, as ex- 
plained in Sec. 6.2.4, there is an important cancellation between the two, neither of 
them can be safely neglected. Without entering here any of the technical details, it 
is difficult to solve the DMFT problem including exactly both these intra-site cou- 
plings, and some approximation must be made. Very detailed DMFT work has been 
done including a realistic representation of the electron-vibrational JT coupling for 
the fuUerides. In particular, Refs. [Ill, 207, 208] study a d = 3 Hubbard- Holstein 
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Fig. 6.16. Superconducting Tc as a function of U for the t®H and a® A models for 
the half filled d = 3 model, treated within the DMFT. The parameters are A = 0.6 
and u)/W = 0.25. The figure illustrates the important difference between Hg and 
Ag phonons. (Prom Ref. [207].) 

model including a single JT Hg phonon mode (5 oscillators rather than 40). The role 
of different lattice symmetry is investigated in Ref. [208]. These approaches employ 
very accurate and extensive quantum Monte Carlo calculations in a finite cluster, 
including the vibronic coupling with realistic JT strength A and a realistic band 
structure. (An earlier calculation [185] even investigated the role of disorder in these 
systems) . Reference [207] addresses mainly superconductivity in the d = 3 Hubbard- 
Holstein-JT model. As illustrated in Fig. 6.16, it was found (like in the d = 2 model) 
that the JT coupling to degenerate phonons produces a pairing much more com- 
patible with strong correlations (large U) than coupling to a single non-JT Ag 
mode. In non-JT Ag electron-phonon coupling superconductivity disappears when 
U„ = U +J7^~^ > 0, while in JT Hg coupling superconductivity survives up to largely 
repulsive Lin > 0. On the other hand, for small U, the local pairing becomes gener- 
ally less efficient, because charge fiuctuations induced by electron hopping disrupt 
the JT ground state(s) into uncorrelated electron pairs. Hence the superconducting 
order parameter is always depressed by U at small U. As U is increased, coherent 
electron hopping is gradually suppressed and the local pair formation becomes more 
important. For realistically large JT coupling A — ^N{eF)^i, glhuik 1 [16, 93] 
d = n = 3 and increasing U the metal phase below Uc is always superconducting, 
with Tc a monotonically decreasing function of U. In other words, electron-electron 
repulsion always disfavors superconductivity for realistically large JT coupling [207] . 

However, this approximation is unrealistic as it neglects Hund's rule exchange, 
and the large cancellation of JT energetics that it implies. Inclusion of exchange 
is difficult, because it leads to a fermion sign problem in the Monte Carlo DMFT 
impurity solver. If on the other hand the impurity is solved by, e.g., Lanczos diagonal- 
ization, exchange can be treated but the phonon ladder gives rise to too many states 
per each site. An approximate way out is available very close to the Mott transition. 
Here the quasiparticle bandwidth zW tends continuously to zero, and sufficiently 
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close to Uc it falls below the typical vibrational frequency tiu). In this regime, the 
frequency dependence of the retarded JT coupling (assumed to be weak) can be inte- 
grated out, resulting simply in an additional but inverted Hund's rule like exchange 
term. As discussed in Sec. 6.2.4, for d = 3 that term is J"'^ = —j^^dk^k- In 
this limit - close enough to the Mott transition, JT coupling not too strong, relar 
tively high vibration frequencies - it is therefore possible to treat both Hund's rule 
exchange and JT coupling, by simply replacing J with JcB — J + J"^"^. This is the 
route generally followed in recent work by M. Capone et al. [146, 186, 190, 191] and 
by M. Granath and S. Ostlund [209]. 

In the tiu molecular level of fullerene J ~ 50 meV, and the DFT-LDA JT 
couplings of Ref. [70] yield J""^ = —57 meV, while the couplings obtained from 
photocmission [93] are consistent with J°^^ = —127 meV: with both estimates the 
net resulting Jeff is negative. The effective Hund's rule of CgJ" is then inverted, 
the low-spin states representing the ground state, as dictated by JT coupling. It 
is important to underline again here that the massive cancellation finally leaves 
fullerides with a weak on-site spin pairing interaction, replacing the original stronger 
couplings, both JT and Hund. The consequences of this weak residual pairing are, as 
it turns out, qualitatively different from either the strong singlet pairing required by 
JT alone, or by the strong triplet pairing required by Hund's J. Similarly, in the half 
filled Cg molecular level of a metal phthalocyanine (ci = n = 2) it is predicted Joff = 
—60 meV [205], again inverted, but again small as a result of a strong cancellation. 

Let us consider in detail the consequences of this cancellation and of the final 
Hund's rule inversion. The first consequence is that the nature of the Mott insulator 
itself is affected. Since the isolated molecular ion ground state has low spin (a spin 
singlet for even n and a spin doublet for odd n), so will the Mott-insulator sites. The 
insulating state is more complicated and intriguing than either the one-band Mott 
insulator, or the simple orbitally-ordered cooperative JT distorted state. It is rather 
a Mott Jahn- Teller insulator. Let us try a description of that state, in the easy case 
of n = 2 electrons/site. Assume initially zero hopping between sites. On each site, 
the JT phonons are characterized by a 2-dimensional (pseudo)rotor in the trough. 
Subject to its own quantum fluctuations, this free rotor favors no special direction 
and will have its (say) L = ground state separated by a gap 5 from its L = 1 
first excited state. Turning now on a hopping oc W between sites, this will cause an 
intersite interaction between rotors of order W'^ /U). If that is strong enough, the 
rotors will freeze statically in an orbitally ordered state, similar to a ferroelectric, 
which is on fact a cooperative JT distortion [76, 102, 103, 210]. If instead the on-site 
quantum fluctuations make 5 large enough, then the static cooperative order will 
quantum mechanically melt. This quantum melted state, the rotor analogue of a 
spin singlet, is the Mott-JT state [211]. This state bears some resemblance to the 
quantum paraelectric state [212]. It is strongly nonadiabatic in nature, which reflects 
in a strongly entangled admixture of electronic and phonon states: all electronic 
spectroscopies should observe radically renormalized and vibronically broadened 
"bands", even at low temperature. Other molecular and lattice properties of this 
state still need to be worked out and investigated. This remains a theoretical task 
for the future. 

Further characterization of the Mott-JT state is needed, both theoretically and 
experimentally, especially when it comes to its spectroscopical properties. The main 

examples at our disposal are NHaKsCeo (orthorhombic) J194] and (NH3)6Li3C60 
(cubic) [12] for n = 3 and a, S = 1/2 on-site ground state, and prominently Rb4C6o, 
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Fig. 6.17. Superconducting gap in a half-filled d = 2-band model, as a function 
of the on-site repulsion for several (anti-Hund) couplings J = Jeff = J + J''~^ 
(where J°~^ accounts antiadiabatically for the electron-phonon interaction if J°~^ — 
— A/[2A''(eF)] = —2g^hLu, and gTiuj represents the coupling strength of the e ® E 
JT coupling at each site; the difference with the relation of Table 6.2 is due to 
d = 2 rather than d = 3). Increasing repulsion spoils superconductivity at strong 
electron-phonon coupling. Superconductivity is instead strongly enhanced close to 
the Mott transition at weak coupling. The inset reports the weak-coupling regime 
on an expanded scale, showing a much smaller gap at small U compared to U ~ Uc- 
(Adapted from Ref. [146].) 

K4C60 for n = 4, 5 = singlet ground state. Concerning the n = 3 Mott insulators, 
the main evidence so far is their clear characterization as spin-1/2 antiferromagnets. 
Since there are 3 electrons/fuUerene, that can be explained only by a JT- dominated 
state. That state could be at this stage either a statically distorted, collective JT 
insulator, or a quantum melted Mott-JT state. The available crystal and magnetic 
structure of orthorhombic NH3K3C60 [168, 194] displays a very rich interplay of 
magnetic and orbital order, which is only well analysed but is apparently compati- 
ble with either possibilities. It should be noted in addition that merohedral disorder 
in the relative angular orientation of the Ceo molecules will frustrate static collec- 
tive order, additionally favoring the melted state. Indeed only in the well ordered, 
"discrete" salts well defined collective distortions have been observed so far [15]. 
An orbitally-ordered static collective JT state [76, 102, 103, 210] has never been 
reported in the alkali fuUerides [213]. On the other hand, in a band-degenerate con- 
text, the possibility of dynamical orbital order should always be considered [211]. 
The question whether one or more of the existing Mott insulating fuUerides could 
be characterized as a dynamical Mott Jahn Teller insulator stands at this stage as 
an exciting experimental challenge. 

After this illustration of the important effects of on site interactions on the mott 
insulating fuUerides, we can now move on to the metallic phase, obtained at smaller U 
values. As we will see, the metallic phase below Uc is even more affected, and in fact it 
often becomes superconducting. This is not in itself totally surprising, since we have 
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Fig. 6.18. Piiase diagram od the d = 2 model as a function of U/W and doping 
5 — n — 2 at Jeff — —0.05 W . The thick vertical line marks the singlet Mott insulator. 
For U = 0.92 W , the inset shows the superconducting gap divided by a factor 3 and 
the Drude weight D (i.e. the strength of the superfluid peak normalized to the 
noninteracting value) as functions of doping. (From Ref. [146].) 



after all a small but definitely attractive Jeff reflecting a locally dominant electron- 
vibration coupling (even if orders of magnitude smaller than the repulsive U\). What 
is more surprising is the behavior of the superconducting gap as a function of U, see 
Fig. 6.17. So long as ( — Jeff) is large - as in the uncompensated JT case described 
earlier - the gap decreases monotonically with increasing U. But when (—Jeff) is 
small - reflecting a strong cancellation by Hund's rule - superconductivity only 
survives near U = and near U = Uc- Superconductivity near the Mott transition 
was designated "strongly correlated superconductivity" (SCS). In SCS, the gap near 
Uc can increase for increasing U, and is many orders of magnitude larger than even 
the repulsion-free BCS-like superconducting gap near [/ = 0. Moreover, close to the 
Mott insulator, the SCS gap magnitude - and thus presumably the superconducting 
Tc - has a characteristic bell-shaped behavior. 

The non-superconducting, metallic phase obtained by sligthly doping the Mott- 
JT insulator at 17 > (7c away from stoichiometry is also strongly unconventional, in 
that it has a pseudogap. Even if doping away from stoichiometry is not straightfor- 
ward in fuUerides, this aspect of the theoretical phase diagram is quite illuminating. 
The presence of the pseudogap is related to Jeff < and thus to the slightly pre- 
vailing JT effect, but is not at all bipolaronic, since as illustrated in Ref. [191] the 
electron occupancy of each site in this state is strongly pinned to the mean value 
n, the bipolaronic disproportionation into (n — 1) and (n + 1) pairs completely sup- 
pressed and in fact reversed by the large Hubbard U. Upon further doping of the 
pseudogap metal, the model exhibits SCS, the maximum gap attained at some op- 
timal doping, then dropping and disappearing in a final overdoped phase, as shown 
in the overall T = phase diagram of the d = 2 n ~ 2 model reported in Fig. 6.18. 

As one can see, many of these features are very reminiscent of high- Tc cuprates, 
whereas the model does not describe cuprates at all. Here we have on-site orbital 
degeneracy, Hund's rule exchange, and JT coupling, all elements that are absent or 
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otherwise quite different in the cuprates. The strongly correlated superconducting 

order parameter in the model is, as in actual fulloridos, s-wavo and not d-wave as in 
the cuprates. We suggest that our models and at least some of the fuUerides which 
they describe, are new members of a wider family of strongly correlated supercon- 
ductors, which includes the cuprates, and all potentially high- Tc materials. 

In this light, it is important to extract and to expose our overall qualitative 
understanding out of the model calculations. The occurrence of superconductivity 
in our orbitally degenerate molecular conductors, meant to model the fuUerides, 
can be qualitatively rationalized by discussing what happens at integer filling as U 
is raised to approach the Mott insulating transition at Uc- First, the quasiparticle 
bandwidth narrows indefinitely, from W to zW with 2 — > as ?7 — > ?7c. The rest of 
the single-particle spectral weight moves to the incoherent Hubbard bands, far away 
from the Fermi level. Second, since all charge fluctuations are gradually frozen out 
when approaching the Mott insulator, so are all the charge-related repulsions (and for 
that matter, attractions too) between the quasiparticles. For example, U is effectively 
reduced to ~ zU , and very near the Mott transition the propagating quasiparticles 
tend to infinitely massive interacting fermions, with a decreasing absolute repulsion. 
Third, thanks to the orbital degeneracy the spin fluctuations are not frozen out, and 
so attractions that act in the spin channel are not renormalized away. Therefore, e.g. 
the effective Hund-rule exchange J or Jeff retains its bare value even at the Mott 
transition, and that can easily overcome the weak repulsion zU . 

These three elements make the effective quasiparticle Hamiltonian very similar 
to an attractive Hubbard model, with a bandwidth zW and an on-site attraction 
Jeff < 0. The ground state is a weak-coupling BCS-like superconductor, at least 
sufficiently below the Mott transition, so long as zW ^ \Jca\- As U grows, the 
superconductivity changes from weak to strong coupling, the maximum gap achieved 
when zW ^ | Jeff |. This point is equivalent to the maximum of the Nozieres-Schmitt- 
Rink [214, 215] curve of Tc versus A in purely electron-phonon superconductors and, 
we believe, also to the optimal doping point in the cuprates. 

If assumed to hold quantitatively, the qualitative mapping onto the attractive 
Hubbard model has the additional virtue of even predicting the maximum 3D crit- 
ical temperature of our type of model at optimal doping, by reading it off the 
attractive Hubbard model studies. That gives approximately fcBT"™^ ~ 0.07|Jefi| 
for our models. In trivalent fuUerides the value of | Jeff| can be extracted by equat- 
ing the observed spin gap ~ 140 meV [127, 128] to 5 | Jeff | (sec Table 6.4), yielding 
jmax ^ Q Qrj I j^giy/jg ~ 23 K, which is surprisingly, though probably unfairly, accu- 
rate. 

One remaining question is: what is the relationship of this "island" of SCS near 
the Mott transition, to the more standard BCS-like phonon driven superconductivity 
one finds in the same model when there is no electron-electron repulsion, i.e. near 
[7 = 0? The answer to this question is very instructive, and is obtained by solving the 
model for a grid of on-site attractive coupling values - negative Je« values [146] . As 
shown by the plot of zero temperature superconducting gaps in Fig. 6.17 if ] Jeff | is 
large, then the pairing gap is in fact maximum at [/ = 0, and decreases monotonically 
for increasing U , vanishing just at the Mott transition. This scenario is, we believe, 
equivalent to that presented by Han et al. for d = n = 2 [216] and d = n = 3 
[111, 207, 208]. In particular, the monotonic decrease of Tc as U approaches Uc in 
Fig. 6.19 (squares) is equivalent to the large-j Jefij monotonic gap of Fig. 6.17 in a 
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Fig. 6.19. (a) Superconducting transition temperature Te vs Coulomb repul- 
sion Z7 in a d = 2-band model with JT electron-phonon coupling, where the cou- 
pling strength ghuj is gauged by J = J""^ = —2g^hui. Weak JT phonon coupling 
{J/W = —0.1) produces superconductivity with a crossover from the conventional 
superconductivity (small U) to the local pairing regime at large U (a Mott-JT in- 
sulator), (b) Quasipaxticle renormalization factor z. (Adapted from Ref. [216].) 

model including explicitly the JT phonons. Similarly, the surge in Tc in Fig. 6.19 
(circles) is essentially identical to that found earlier for small attractive |Jeff|. 

Things do change qualitatively when the effective attraction decreases, as the 
dominance of JT is more and more eroded and cancelled by Hund's rule exchange. 
For decreasing IJcffj, superconductivity remains strongest at [/ = and near the 
Mott transition, but it gradually weakens in between. Eventually for very small 
I Jeff I one arrives at two separate superconducting islands, one with a tiny BCS gap 
near U = 0, and a second one, near the Mott transition at Uc, with a gap orders of 
magnitude larger than BCS. This the SCS island. At the lower edge of this island 
(equivalent to the "overdoped" regime of the cuprates), the pairing is again BCS 
like: but the gap and presumably Tc rises as the electron-electron repulsion increases, 
rather than the opposite. We believe that this regime could describe most of the 
trivalent fullerides, and their gigantic Tc and susceptibility increase with volume 
expansion. In fact the whole curve of Tc versus increasing volume as presented, e.g., 
by Durand et al. [12] (Fig. 6.20) could in our view be roughly identified with the 
theoretical bell-shaped curve of gap versus increasing U of Fig. 6.17. 

A corollary is that the highest Tc fullerides should be extremely strong-coupling 
materials, comparable to cuprates. This is in agreement with the extremely large 
Hc2 values observed experimentally as well as with the existence of an irreversibility 
line in the H — T phase diagram [217]. Both features are presently unexplained and 
both are similar to cuprates. 
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Fig. 6.20. A schematic pliase diagram of CgQ fullerides. The superconducting 
transition temperature Tc is shown as a function of the ratio U/W divided by the 
critical value Uc/W required to produce electron localization - Uc/W is estimated to 
be 2.3 for the f.c.c. phases (square and diamond symbols), 1.3 for b.c.c. CsaCeo and 
Lis (NH3)6C6o (circles) and 1.1 for the orthorhombic NH3 AaCeo- The value of U 
used to calculate U /W was 0.8 eV. The Mott-insulator region is shaded: uncertainties 
in the precise numerical estimates are emphasised by the graded shading of the 
crossover between metallic/superconducting and localized. W is derived for the f.c.c. 
phases from Ref. [2], for the orthorhombic NH3 K3C60 from Ref. [140] and for b.c.c. 
phases by DFT calculation [12]. The dramatic effect of lifted orbital degeneracy 
on the location of orthorhombic NH3 K3C60 is illustrated by comparison with an 
hypothetic "cubic" NH3K3C60) retaining the tiu orbital degeneracy and located 
well within the metallic regime. (From Ref. [12].) 



The compounds in the descending Tc branch Fig. 6.20, apparently very close to 
the Mott transition, are the fuUeride analog of the "underdoped" cuprates. They 
should be especially interesting, and their normal state should for example exhibit 
strong pseudogap features. 

So far the discussion assumed prevalence of JT over Hund's rule exchange, JcB < 
0. Also a hypothetical case of positive Jeff (regular Hund's rule) would not be without 
interest. On-site high-spin states would prevail in this case. Again, as compared to 
J = 0, the Mott transition shifts to lower Uc [216]. This occurs for two reasons: a 
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reduction of the relevant pair energy Un proportional to J, and the fact that splitting 

of different Hund's multiplcts opposes the metallic state. In the DMFT calculation 
of Ref. [216], the possibility of triplet superconductivity was investigated, and indeed 
it was found that before the metallic state turns into a magnetic Mott insulator, 
a superconducting phase emerges, with triplet pairing, but without p wave, thanks 
to the orbital degeneracy. If for example CgJ was really high spin in (AsF6)2C6o 
or (SbF6)2C60, among the metallic phases obtainable under pressure one could find 
one such triplet SCS. 

In practice a small negative anti-Hund Jefi represents better the net balance of 
the positive Coulomb term and a prevailing electron-phonon attraction (assumed 
in the antiadiabatic limit) in C^^ [105, 106, 112]. Figure 6.17 shows that a strong 
reduction of Uc occurs under the action of an even small negative Jeff (as expected 
for any local interaction). In addition, as discussed in Sec. 6.2.4, a negative Jeff 
favors local low-spin states, which are spin singlets whenever the local occupancy 
is even. Of course, in the solid, band effects tend to favor local charge fluctuations 
and spin admixture, but as U approaches Uc, the bandwidth reduces effectively to 
the width zW oi the Kondo peak. When z W becomes comparable to | J| , the local 
exchange term promotes locally paired singlets in a sort of anti-Migdal regime [176, 
177, 178, 179, 180, 181] close to the transition, where the average phonon frequencies 
exceed the effective bandwidth. This regime favors superconductivity, and DMFT 
suggests that a high-Tc (fcBlc of order | Jeff |) superconducting phase should intrude 
between the metal and the low-spin Mott insulator [146, 186, 190, 191]. As shown in 
Fig. 6.17, for very small [Jeff] ^ W (representative of the cancellation characteristic 
of Cg,7 ) , the very weak superconductivity of the uncorrelated regime C/ ~ becomes 
dominating close to the Mott transition. An alternative possibility to the intruding 
superconducting phase is a discontinuous, first-order transition, jumping directly 
from the metal to the insulator, and in fact this is a serious eventuality for the A4C60 
fullerides, where no superconductivity is observed. The pressure-induced reversion 
from insulator to metal observed in tetragonal Rb4C6o [126] is most likely the n = 4 
counterpart of the n = 3 transition of the ammoniated compounds. 



6.5 Conclusions 

We have discussed some results and qualitative theoretical ideas on the modeling 
of metal-insulator transition in fuUerene compounds. They have been described as 
Hubbard models with electron hopping with a bandwidth W between orbitally de- 
generate fuUcrene sites, each of them supporting a Coulomb repulsion U, a Hund's 
rule exchange J, and a JT distortion with energy gain E''~^ . The interplay of these 
parameters gives rise to a rich phase diagram, comprising Mott insulators, metals, 
and superconductors. For electrons in fullerene, the JT strength marginally prevails 
over exchange, a fact which has important implications for all phases. The Mott 
insulators can in reality become Mott-JT insulators, whose local configuration is 
low spin, with an electron-vibrational entanglement as intrinsic as that of a su- 
perconductor, and possibly displaying other unexplored features. The metal phases 
obtained by doping the even-n Mott insulators can possess a pseudogap, namely a 
depression of the density of states at the Fermi level, and a very reduced suscepti- 
bility, constituting a kind of nonadiabatic semimetals. The superconductor phases 
near the Mott insulators arise through pairing of quasiparticles, these constituting 
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a thin conducting web floating in a Mott-insulating sea. These model superconduct- 
ing phases, ahhough ultimately caused by JT effect (electron-phonon) and s-wave 
in character, are shown to share many features with those of high- Tc cuprates. This 
possibility calls for further experimental and theoretical work. 

Further work should also be devoted to investigate the detailed nature of the 
Mott phase transition. The alternative possibilities of first-order and second-order 
transition, should be distinguised in the n = 3 and n — 4 classes of fuUerides. For 
n = 3 compounds, the simultaneous disappearence of supercondtucting order and 
appearence of magnetic order seem theoretically incompatible with a continuous 
transition. 

To investigate possible spin/orbital antiferro correlations of the kind observed in 
NH3 K3C60 [194], one should study a detailed three-dimensional model, including the 
highly directional hopping matrix elements between different Ceo sites (and possibly 
merohedral disorder). Some of these aspects could be addressed by standard mean- 
field techniques (of the type applied in a different context [218]), or quantum Monte 
Carlo cluster calculations which unfortunately suffer by significant finite-size effects. 
The DMFT method, based on an infinite lattice is free from these limitations, but 
the same aspects are inaccessible to a single-impurity DMFT model: cluster DMFT 
methods [183, 184] could lead to some progress in this direction. 
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